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Hydrodynamic transport coefficients may be evaluated from first principles 

ON 



in a weakly coupled scalar field theory at arbitrary temperature. In a theory 
with cubic and quartic interactions, the infinite class of diagrams which con- 



tribute to the leading weak coupling behavior are identified and summed. 



The resulting expression may be reduced to a single linear integral equa- 
tion, which is shown to be identical to the corresponding result obtained 
from a linearized Boltzmann equation describing effective thermal excita- 
tions with temperature dependent masses and scattering amplitudes. The 
effective Boltzmann equation is valid even at very high temperature where 
the thermal lifetime and mean free path are short compared to the Compton 
wavelength of the fundamental particles. Numerical results for the shear and 
the bulk viscosities are presented. 
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I. INTRODUCTION 



Linear response theory provides a framework for calculating transport coefficients start- 
ing from first principles in a finite temperature quantum field theory. As reviewed below, 
the resulting "Kubo" formulae express the hydrodynamic transport coefficients in terms of 
the zero momentum, small frequency limit of stress tensor-stress tensor correlation func- 
tions |TJ. One- loop calculations of transport coefficients using these Kubo formulae in a 
relativistic scalar 4 theory have appeared previously However, those calculations 

are wrong even in the weak coupling limit; they fail to include an infinite class of diagrams 
which contribute at the same order as the one-loop diagram. These multi-loop diagrams 
are not suppressed because powers of the single particle thermal lifetime compensate the 
explicit coupling constants provided by the interaction vertices.^] 

In this paper, all diagrams which make leading order contribution to the viscosities 
in a weakly coupled relativistic scalar field theory with cubic and quartic interactions are 
identified. The diagrammatic rules needed to calculate the required finite temperature 
spectral densities of composite operator correlation functions were derived in a previous 
paper []5| (and are summarized below). The dominant diagrams are identified by counting 
the powers of the coupling constants which result from a given diagram, including those 
generated by near "on-shell" singularities which are cut-off by the single particle thermal 
lifetime. 

For the calculation of the shear viscosity, certain cut "ladder" diagrams, corresponding 
to the contribution of elastic scatterings only, are found to make the leading order contribu- 
tions. The geometric series of cut ladder diagrams is then summed by introducing a set of 
effective vertices which satisfy coupled linear integral equations. The resulting expression 
is then shown to reduce to a single integral equation, which is solved numerically. 

For the calculation of the bulk viscosity, in addition to the leading order ladder dia- 



Similar phenomena occur in the calculation of transport coefficients in non-relativistic fluids 
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grams, contributions from the next order diagrams containing inelastic scattering processes 
must also be summed. In general, the bulk viscosity is proportional to the relaxation time 
of the processes which restore equilibrium when the volume of a system changes || . For a 
system of a single component real scalar field, such processes involve inelastic scatterings 
which change the number of particles. Hence, diagrams corresponding to such processes 
must be included. 

Boltzmann equations based on kinetic theory have traditionally been used to calculate 
transport properties of dilute weakly interacting systems. However, the validity of kinetic 
theory is restricted by the condition that the mean free path of the particles must be much 
larger than any other microscopic length scale. In particular, the mean free path must 
be large compared to the Compton wavelength of the underlying particle in order for the 
classical picture of particle propagation to be valid. A Boltzmann equation describing the 
fundamental particles cannot be justified when this condition fails to hold. Such is the 
case at extremely high temperature, where the mean free path scales as 1/T. 

No such limitation exists when starting from the fundamental quantum field theory. 
Nevertheless, it will be shown that the correct transport coefficients, in a weakly coupled 
theory, may also be obtained by starting from a Boltzmann equation describing effective 
single particle thermal excitations with temperature dependent masses and scattering am- 
plitudes.0 This equivalence holds even for asymptotically large temperatures where both 
the thermal lifetime and mean free path are tiny compared to the zero temperature Comp- 
ton wavelength. Hence, in a weakly coupled theory, although a kinetic theory description 
in terms of fundamental particles is only valid at low temperatures, a kinetic theory de- 
scription of effective thermal excitations remains valid at all temperatures.^ 



2 The Boltzmann equation derived by Calzetta and Hu j?]] via relativistic Wigner function is also expressed 
in terms of the thermal mass. 

3 In a strongly coupled theory, the mean free path can be comparable to the scattering time, or other 
microscopic scales, and no kinetic description is justified. 
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The effective kinetic theory result presented in this paper is valid through all tem- 
peratures in the weak coupling limit. At low temperatures, T <C m phys , where m p h ys is 
the physical mass of the underlying particles at zero temperature, most particles are non- 
relativistic. Hence, the effective theory reduces to non-relativistic kinetic theory at low 
temperatures. If the temperature is in the range m phys < T <C m phys / y/\, where A is the 
quartic coupling constant, most particles are relativistic, but the thermal corrections to 
the mass and the scattering amplitude are negligible. Consequently, the viscosities at these 
temperatures can be calculated by a kinetic theory of relativistic particles with temperature 
independent mass and scattering amplitudes. 

The most interesting temperatures are those where T = 0{m v ,i iys /\/~X). At these tem- 
peratures, the thermal correction to the mass is comparable to the zero temperature mass. 
For weak coupling, this temperature is also large enough that most excitations are highly 
relativistic. One might expect that the thermal correction to the mass would then be ir- 
relevant. This is true for some physical quantities which are insensitive to soft momentum 
contributions, such as the shear viscosity. However, other quantities, such as the bulk vis- 
cosity, are sensitive to soft momenta. For such quantities, including the thermal correction 
to the mass and the scattering amplitude will be shown to be essential. 

At very high temperature, T > m p h ys /A, all mass scales, including the cubic coupling 
constant, other than the temperature are completely negligible, and consequently the the- 
ory reduces to the massless scalar theory with only a quartic interaction. 

Through out this paper, we work with the Lagrangian 

-£ = \<K-% - v 2 + ml)<P + f 3 + ^ • (1.1) 

It is assumed that A 1 and g 2 = 0(Xm 2 ) ), so that the theory is always weakly coupled. 
For simplicity, we also take m„ > 0. Note that at tree level, m can be regarded as 
the physical mass m phys . Portions of the analysis will begin by assuming pure quartic 
interactions, after which the additional contribution arising from cubic interactions will be 
considered. The remainder of the paper is organized as follows. A brief review of various 
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background material is presented in section ||. This material includes the definition of 
transport coefficients, basic linear response theory, diagrammatic "cutting" rules for the 
evaluation of spectral densities, and a summary of the behavior of self-energies at high 
temperature. Section [TTI| deals with the problem of identifying the leading order diagrams. 
By counting powers of coupling constants, including those from near on-shell singularities, 
ladder diagrams are identified as the leading order diagrams. The summation of these 
diagrams is discussed in section |TV|. Section [V| contains a brief review of the computations 
of viscosities starting from the Boltzmann equation, and then discusses the relation between 
the resulting formulae and those in section |TV|. Using the results of the previous sections, 



the final calculation of viscosities is discussed in section [VT|, and numerical results presented 
in section [VII] . 

Several appendices contain technical details. Explicit forms of the imaginary-time and 
real-time propagators used in the main body of the paper are summarized in appendix [A]. 
Appendices (B| and |C| present explicit forms of the "ladder" kernels discussed in section [TV[ 
In appendix [B|, the first order correction to the equilibrium stress-energy tensor needed in 
sections [TV] and [VT] are calculated. Appendix [E] discusses the soft momentum and collinear 
contributions to finite temperature cut diagrams, and shows that they does not upset the 
estimates used in section |HI[ Appendix [F] contains technical details of summing up the 
"chain" diagrams appearing in section [I I [~D| . 



II. BACKGROUND MATERIAL 

A. Hydro dynamic transport coefficients 

In a single component real scalar field theory, the only locally conserved quantities are 
energy and momentum. The transport coefficients associated with energy and momentum 
flow, known as the shear and bulk viscosities, may be defined by the following constitutive 
relation, 
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Fa) * (V,^ ) + VjQi) - ¥aV l (T?)) - ^-£— Sij V l (lf) + 6v(V) , (2.1) 



valid when the length scale of energy and momentum fluctuations is much longer than the 
mean free path. Here, is the spatial part of the stress-energy tensor T^, e = T 00 is the 
energy density, V = |T? is the pressure, and r\ and ( are the shear and the bulk viscosities, 
respectively. Also, (■ ■ ■) denotes the expectation in a non-equilibrium thermal ensemble 
describing a system slightly out of equilibrium. Since there are no additional conserved 
charges, thermal conductivity is not an independent transport coefficient .[] 

The above constitutive relation, combined with the exact conservation equation 

d,T^(x) = , (2.2) 

constitute linearized hydrodynamic equations for a relativistic fluid. With the help of the 
equilibrium thermodynamic relation 

d(s) C(l ~ V ° ' [2J} 

where v s is the speed of sound, the linearized hydrodynamic equations can be reduced to 
two equations for the transverse part of the momentum density (tt t ) and the pressure (V), 

(d t - DV 2 ) (it t (x)) = , (2.4a) 

and 

(d 2 t - v 2 s V 2 - TV%) (V(x)) = . (2.4b) 
Here the diffusion constant D is proportional to the shear viscosity 

D = v /(e+V) eq , (2.5) 

and the sound attenuation constant T equals a linear combination of the viscosities 

T = (lr ]+ C)/(e+V) eq . (2.6) 



4 Calculation of the thermal conductivity in a scalar \<fi 4 theory in Ref. S is in this sense misleading. 
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Using the basic linear response result, one may express the viscosities in terms of the 
stress tensor-stress tensor correlation functions []I]J§. One finds the "Kubo" formulae, 

71 = ^ lim n ^^(^ <j) ' ( 2 - 7a ) 
20 "J-»0 q^O 

Here a^ n (u, q) is the Fourier transformed traceless stress-stress Wightman function, 

a^(u,q) = J d 3 xdte~^ + ^ (n lm (t,*)n lm (0)) eq , (2.8) 

where 

7r lrn (x) = T lm (x) - \5i m T&x) (2.9) 
is the traceless stress tensor. Similarly, 

app(u,q) = J d z *dte- i(iM (V(t, x)P(0)) cq , (2.10) 

where 

P(t,x) = V(t,x)-v 2 e(t,x) = \Tl{x) - v 2 T m {x) , (2.11) 

is a linear combination of pressure and energy density. The constant v 2 in this combina- 
tion is arbitrary; due to energy-momentum conservation, Wightman functions involving 
the energy density vanish (at non-zero frequency) in the zero spatial momentum limit. 
However, as will be discussed in section [TV], v will eventually be chosen to equal the speed 
of sound. This will be necessary in order to make the final integral equation for the trans- 
port coefficient well defined. Note that, whereas the approximate constitutive relation 
( p.l|) involves a non-equilibrium thermal expectation, the Kubo formulae (|2.7|) express the 
transport coefficients solely in terms of equilibrium expectation values. 



B. Qualitative behavior of the viscosities 



In general, a transport coefficient is roughly proportional to the mean free path, or 
equivalently the relaxation time, of the processes responsible for the particular transport 



This behavior is most easily seen in a diffusion constant (or in the shear viscosity). 
Consider a system with a conserved charge. In such a system, the diffusion of a charge 
density fluctuation may be modeled by a random walk || . The rate of the diffusion then 
depends on two parameters; the step size (the mean free path) and the number of steps 
per time (the mean speed). A longer step size or a larger number of steps per time implies 
faster diffusion of the excess charge, i.e., a larger diffusion constant D. Since the diffusion 
constant has the dimension of a length, one finds 

D~l iTec v. (2.12) 

Recall that the diffusion constant in Eq. ( |2.5| ) is given by D = r// (e+V) eq . Applying the 
above estimate of D yields 

V ~ ^frcc^+^eq • (2.13) 

Given the scattering cross-section a and the density of the particles n, the mean free 
path can be estimated as Zf ree ~ 1/na. Consider a weakly coupled scalar A0 4 theory. The 
lowest order scattering cross section in the A0 4 theory is a ~ X 2 /s where s is the square of 
the center of mass energy. At high temperature, T 3> m p h y s, the only relevant mass scale 
is the temperature. (Here m p h ys denotes the physical mass.) Hence, a ~ A 2 /T 2 , n ~ T 3 , 
and 

r/~T 3 /A 2 (T>m phys ). (2.14) 

At low temperature, T <C rn piiys , a ~ A 2 /mp hys and v ~ (^/^phys) 1 ^ 2 - At these temper- 
atures, the energy density (e) eq ~ rn phys n, dominates over the pressure. Canceling two 
density factors in Zf ree and (e) eq , the shear viscosity can be estimated as 

V ~ m 3 hys (T/m phys ) 1 /2 /A 2 (T <<; mphys) (215) 

Note that the shear viscosity is not analytic in the weak coupling constant. This may be 
taken as an indication that the first few terms in the usual Feynman diagram expansion 
cannot produce the correct value of the leading order shear viscosity. 
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For the bulk viscosity, the situation is more complicated than the simple picture given 
above. The bulk viscosity does not have an interpretation as a diffusion constant. Hence, 
the random walk model cannot be directly applied. The bulk viscosity is still proportional 
to the mean free time (inverse transition rate) of a scattering process since the viscosities 
govern relaxation of a system towards equilibrium. However, the factors multiplying the 
mean free time cannot simply be (e+V) eq since the bulk viscosity ( vanishes in a scale 
invariant system ||. 

To understand this, consider a slow uniform expansion of the volume of a system. In 



such an expansion, there can be no shear flow ||10|| . Hence, the relaxation of disturbances 
caused by the expansion depends only on the bulk viscosity. For scale invariant systems, 
the restoration of local equilibrium does not require any relaxation process. A suitable 
scaling of the temperature alone can maintain local equilibrium at all times. Hence, for 
such systems, including the non-relativistic monatomic ideal gas0 and the ideal gas of 
massless particles, the bulk viscosity vanishes since the relaxation time vanishes. 

When the system is not scale invariant, the bulk viscosity must be proportional to a 
measure of the violation of scale invariance, or the mass m p h ys . In section [TV], the formula 
for the leading order bulk viscosity at high temperature is shown to be 

C~^hys T free- (2.16) 

The mean free time Tf ree here is given by the inverse of the transition rate per particle, 

r free = n/(dW/dVdt) , (2.17) 

where dW/dVdt the transition rate per volume, corresponding to the relaxation of the 
uniformly expanding system. In a number non- conserving system with broken scale invari- 



5 The non-relativistic monatomic ideal gas is not equivalent to the low temperature limit of the single 
component real scalar field theory. The number of particles in the non-relativistic monatomic ideal gas is 
conserved whereas the number of particles in the low temperature limit of the scalar field vanishes as the 
temperature goes to zero. Hence, the low temperature limit of the scalar theory bulk viscosity need not 
vanish. 
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ance (such as a massive scalar theory), the number- changing inelastic scattering processes 
are ultimately responsible for relaxation towards equilibrium. As the system expands, the 
temperature must decrease since the system loses energy pushing the boundary. Decreasing 
energy implies decreasing particle number in a number non- conserving system with broken 
scale invariance. Hence, the relaxation toward equilibrium must involve number- changing 
scatterings. In the g0 3 +A0 4 theory, the lowest order number- changing process is 2-3 scat- 
terings involving 3 cubic vertices or 1 cubic and 1 quartic vertices. In pure A0 4 theory, the 
lowest order number changing process is 2-4 scatterings involving 2 quartic vertices. 

At high temperature, m p h ys <C T < m p h ys / V% most particles have momentum of 0(T). 
However, due to the Bose-Einstein enhancement, the transition rate per volume will be 
shown to be dominated by the 0(m t h) momentum components in the system where m t h is 
the thermal mass containing 0(\f~\T) thermal corrections. In the g0 3 +A0 4 theory, with the 
statistical factors for 5 particles involved in the scattering, the transition rate per volume 
is 0(X 2 g 2 T 5 /m 3 h ) which is C(T 3 /mf h ) larger than the transition rate of the particles with 
0(T) momentum. Hence at temperatures in the range m p h ys < T < rn phys /^/\, the bulk 
viscosity will be 

C ~ m phys m 3 th /\ 2 g 2 T 2 . (2.18) 

When T = 0(m phys /V\), ( = 0(T 3 /V\) which is C(A 3/2 ) smaller than the shear viscosity. 
At the same temperature, the A0 4 theory transition rate per volume is (9(A 3 T 4 ) again due 
to the Bose-Einstein enhancement. In this case, ( = (D(j3m phys / \ 3 ) = 0(T 3 /\) which is 
0(A) smaller than the shear viscosity. 

At very high temperature T m phys /A, all mass scales, including the cubic cou- 
pling constant, other than the temperature are completely negligible, and consequently 
the theory reduces to the massless scalar theory with only a quartic interaction. The 
massless scalar theory is classically scale invariant. However, quantum mechanics breaks 
the scale invariance. The measure of the violation of scale invariance in this case is the 
renormalization group /3-function. Since the transition rate per volume must still be 
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(9(A 3 T 4 ) due to the thermally generated O(VXT) mass, the bulk viscosity in this case 
is £>(T 3 /3(A) 2 /A 3 ) = £>(AT 3 ). 

At low temperature, T <C ^ p h y s, the bulk viscosity is proportional to the mean free 
time and the number density. At these temperatures, the transition rate per volume in the 
g0 3 +A0 4 theory is 0(e~ 3f3mphys ), since the center of mass energy must exceed 3m p h ys for 
a 2-3 process to occur. Since the density n = (9(e~^ mphys ) at low temperature, the bulk 
viscosity is then 0(e l3mphys / A 3 ). The bulk viscosity at low temperature is hence much larger 
than the 0(1/ X 2 ) shear viscosity. In the A0 4 theory, the transition rate is 0(e~ At3mphys ) 
since the center of mass energy in this case must exceed 4m p h ys for a 2-4 process to occur. 
The bulk viscosity is then (9(e 2/3mphys /A 4 ), and again much larger than the shear viscosity. 

C. Linear response theory 

Linear response theory describes the behavior of a many-body system which is slightly 
displaced from equilibrium. First order time dependent perturbation theory implies that 

E 

5<i a (t,x)> = i / dV/* dt' ([A a (t,x),A b (t',x')]) eq F b (t',x') , (2.19) 

where {F b (t, x)} is some set of generalized external forces coupled to the interaction picture 
charge density operators {A b (t,x)} so that 

SH(t) = - J d 3 x F b {t, x) A b (t, x) , (2.20) 

and (• • -) cq denotes an equilibrium thermal expectation. (Summation over the repeated 
index b should be understood.) 

To examine transport properties, it is convenient to consider a relaxation process 
in which the external field is held constant for a long time (allowing the system to re- 
equilibrate in the presence of the external field), and then suddenly switched off, 

F 6 (t,x) = F b (x) e et 6(-t) , (2.21) 
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where e is a positive infinitesimal number. Once the field is switched off, the system 
will relax back towards the original unperturbed equilibrium state. Spatial translational 
invariance implies that Fourier components of the initial values 8(A a (0, x)) are linearly 
related to the Fourier components of Fb(x). Hence, after a Fourier transform in space and 
a Laplace transform in time, Eq. ( 2.19Q turns into an algebraic relation ||, 



5A a (z, k) = hx ab (z, k) Xbc(^ k) - 5 ac }5A c (t=0, k) , (2.22) 

%z 

where 5A a {z ) \s) are Laplace and Fourier transformed deviations from equilibrium values 
S(A a (t, x)), and SA c (t=0,k) are Fourier transformed initial values 5(A c (t=0, x)). Here, 
Xabi z i k) is the retarded correlation function with complex frequency z; it has the spectral 
representation 

M)= rtoe*m, (2.23) 
J Ztt uj — z 

where the spectral density is 

Pab (u,k) = J Ae" w ([i o (t,x),i 6 (0)]) eq . (2.24) 

If AaS are conserved charge densities, then Ward identities can be shown to imply 
that the response functions have hydrodynamic poles (poles in the frequency plain which 
vanish as the spatial momentum goes to zero) |TTJ. In the case of the conserved energy 
and momentum densities, the response functions in Eq. (|2.22j) can be shown to have a pole 
at z = —iD\t 2 when the disturbed charge is the transverse part of the momentum density 
7r T , and poles at z 2 = v 2 \t 2 — iTzk 2 when the disturbed charge is the energy density or the 
longitudinal part of the momentum density. 

Eq. ( p.22j ) solves the initial value problem in terms of the response function. When 
the conserved quantities are energy and momentum densities, the time evolution of the 
initial values can be also described (for low frequency and momentum) by the phenomeno- 
logical hydrodynamic equations Q2.4p . When Fourier transformed in space and Laplace 
transformed in time, Eqs. (|2.4j) yield response functions with exactly the same diffusion 
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and sound poles. Hence, by extracting the diffusion constant D and the sound attenuation 
constant T from the pole positions in the correlation functions, one may derive the Kubo 
formulae ( |2.7| ) for the viscosities. 

The Wightman functions appearing in formulae ( p.7|) for the viscosities are trivially 
related to the corresponding spectral densities: 

Pmr (u>, q) = (1 - e"*") q) (2.25) 

ppp(u, q) = (1 - e"^) app(u, q) . (2.26) 

Hence, the viscosities can equivalently be written as zero frequency derivatives of spectral 
densities, 

1 d 

V = on hm a hm n -*-P™( u ' ^ ' (2.27a) 

20 U3-*0 q^O QijJ 



and 



1 d 

C = - lim lim —p pp (u, q) . (2.27b) 

2 kj^o q^o o»u; 



D. Cutting rules 

In Ref. H, diagrammatic cutting rules for the perturbative calculation of the spectral 
density of an arbitrary two-point correlation function were derived starting from imagi- 
nary time finite temperature perturbation theory. These rules are a generalization of the 
standard zero-temperature Cutkosky rules, to which they reduce as temperature goes to 
zero. 

To calculate the perturbative expansion of a finite temperature spectral density, one 
should draw all cut Feynman diagrams for the two-point correlation function of interest .0 



6 Only half the cut diagrams, those in which the external momentum flows into the shaded region, need to 
be considered if one includes an additional overall factor of (1— e~ q where q° is the external frequency. 
Omitting this factor, the same rules generate the Wightman function instead of the spectral density. 
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FIG. 1. A typical cut diagram in a scalar Xcj) 4 theory. 

All cuts that separate the two external operators are allowed at non-zero temperature. 
Each line corresponds to either a cut or uncut thermal propagator, as described below. 

A typical example of a finite temperature cut diagram is shown in Fig. |]. Note that 
a cut at finite temperature can separate a diagram into multiple connected pieces, some 
of which are disconnected from the external operators. A disconnected piece, such as the 
portion labeled A in Fig. p], cannot contribute in a zero temperature cut diagram because 
of energy momentum conservation. For example, at zero temperature the piece labeled A 
would represent an impossible event of four incoming on-shell physical particles scattering 
and disappearing altogether. However, at finite temperature there exist physical thermal 
excitations in the medium. Thus, the above disconnected piece also represents the elastic 
scattering of a particle off of a thermal excitation already present in the medium. This 
scattering process is clearly possible; the amplitude is proportional to the density of the 
thermal particles (as the form of the cut propagator shown below clearly indicates). 

An uncut line in the unshaded region corresponds to a real-time time-ordered propaga- 
tor (T(0(x)0(O))), an uncut propagator in the shaded region is (T(0(x)0(O)))*, and a cut 
line corresponds to the Wightman function (<p(x) (f)(0)). In momentum space, the uncut 
propagator has the following spectral representation 
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where p(cu, |k|) is the single particle spectral density. The cut propagator is proportional 
to the single particle spectral density 

S(k) = [l+n(k )] p(k) , (2.29) 

where n(k°) is the Bose statistical factor l/(e fc0/3 — 1). In more physical terms, 

dV(k) ee #(A;°) ^ = #(/c°) [l+n(k°)] p(k) -0^ (2.30) 

is the thermal phase space volume available to a final state particle in a scattering process, 
and 

dN(k) ee 0{-k°) S(k) ^ = 6(k°) n(k°) p(k) ^ (2.31) 

is the number of thermal excitations within the 4- momentum range (k^, k^+dk^). 

If the single particle spectral density is approximated by a delta function (to which it 
reduces at zero temperature), i.e., 

p ilcc (k) = S gn(k°)2rr5(k 2 + ml), (2.32) 

then self-energy insertions on any line generate ill-defined products of on-shell delta func- 
tions. Although these on-shell singularities disappear when all cut diagrams are summed, 
it is far more convenient to first resum single particle self-energy insertions. The resummed 
single particle spectral density p(k) will then include the thermal lifetime of single particle 
excitations, which will smear the ^-function peaks and produce a smooth spectral den- 
sity. Henceforth, all single particle propagators will include the thermal self-energy and no 
self-energy insertions will appear explicitly in any cut diagram. 

E. Propagators and self-energies at high temperature 

To analyze near infrared singularities, the explicit forms of single particle propagators 
will be needed. The resummed single particle spectral density can be calculated as 
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(a) (b) (c) 

FIG. 2. One-loop self-energy diagrams in a scalar g^+Xcj) 4 theory. At high temperatures, the diagrams 
(a) and (b) produce a thermal mass squared of order AT 2 . The contribution of the diagram (c) to the 
thermal mass is C(g 2 T/m th ) = 0(A 3/2 T 2 ). 

the discontinuity of the analytically continued imaginary-time propagator across the real 
frequency axis, 



p(k) = -i (G E (k°+ie, k) - G E (k°-ie, k) 



+ 



k 2 + m 2 h + J]{k) k 2 + m 2 h + E(k)* 
= _ iGsW+iG ^) = _^L_, (2 .33) 

where the subscript R indicates the retarded propagator given by the analytic continuation 
of the Euclidean propagator 

G R (k)=G E (k°+ie,k) , (2.34) 

and the subscript A indicates the advanced propagator defined similarly, but with k°—ie 
instead of k°+ie [|l|,[n|]. Here, the thermal mass m t h includes the 0(XT 2 ) one-loop correc- 
tions shown in Fig. [| and is the analytically continued single particle self-energy 

£(Jfe) = Z E (k°+ie, k) = E i? (A;)-zE J (A;) . (2.35) 

The thermal mass squared m 2 h may be defined by the (off-shell) condition S#(0) = 0. 
Also, note that since p(k) is the spectral density of a correlation function of the CPT even 
Hermitian operators 0, p(k) must be an odd function of the frequency k° j5|||. This implies 
that is also an odd function of k°. 

As will be reviewed below, the imaginary part of the self-energy is 0(X 2 ), and so is 
small for weak coupling. Hence, the definition ( p.33|) shows that the spectral density in 
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the weak coupling limit has sharp peaks near k° = ±Ek, where the effective single particle 
energy E k satisfies the dispersion relation 

E 2 k = k 2 + m 2 th + Z R (E k , |k|) . (2.36) 

Near the peaks, the spectral density may be approximated by a combination of two 
Lorentzians, 

**> = k (w=wrn ~ wrwrn) • (2 - 37) 

where r k is the momentum dependent thermal width given by 

r fc = £ 7 (£ fc ,|k|)/2£ fc . (2.38) 

The thermal width r\ is always positive since Hj(k), or equivalently the single particle 
spectral density, must be positive for positive frequencies. This can be easily seen from 
the relation between the spectral density and the Wightman function and the positivity of 
the (Fourier transformed) Wightman function (0(rc)0(O)) [0. Note that altogether p(k) 
has four poles at k° = E^iTk and k° = —Ek±iTk- In terms of the single particle spectral 
density, the cut propagator is 

„,,, . ,,„,, , lN 2\l+n(k°)}J: I (k) 

s(i) = N ' Mi) = ^f ' (2 - 39) 

For the uncut propagator given by Eq. (|2.28|) , the frequency integral can be exactly 
carried out to yield (see appendix |A| for details) 



|P + m t 2 h + £(A;)| 
In the weak coupling limit, this becomes 

1 + n(E k ) . n(E, 



k 2 + m 2 h + E(k) k 2 + m t 2 h + E(fc)* 
k 2 + m 2 h + Z R (k) coth(A; /j/2) 

2 + o P( k ) ■ (2.40a) 



The first term in Eq. (|2.40b ) has poles at k° = ±(Ek—iTk), and the second term has poles 



at k° = ±(Ek+iTk), coinciding with the pole positions of the spectral density p(k). 
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(a) (b) (c) (d) 

FIG. 3. Cut two-loop self-energy diagrams in a scalar g(j) 3 +\(f> 4 theory which produce a thermal width 
of order X 2 T. 

An important point to notice is that even though the statistical factor n(k°) has a pole 
at k° = 0, both the cut propagator S(k) and the uncut propagator G{k) are finite at zero 
frequency since the self-energy £/(&), which is an odd function of k°, vanishes at k° = 0. 
Hence, although numerous factors of the statistical factors may appear in an expression 
for a diagram, one can be sure that there is no pole when loop frequencies approach zero. 

To determine the size of the thermal mass m t h and the thermal width at high 
temperatures, the size of the one-loop (Fig. |2|) and the two-loop (Fig. |3]) self-energies at 
on-shell momenta must be known. At relativistic temperatures, T > m , the first one- loop 
diagram, Fig. Qa, generates an 0(XT 2 ) contribution to the real part of the self-energy. 
Diagram with two cubic interaction vertices, is 0(g 2 T 2 /m 2 h ) which is at most 0(XT 2 ) 
since by assumption g = O(VXm ), and m t h > rn . At high temperature, the real part of 
diagram Fig. |c is 0(g 2 \n(T/m th )) which is at most 0(X 2 T 2 In A). Hence, the diagram |c 
does not contribute to the leading weak coupling behavior of the thermal mass correction]] 
Hence, the thermal mass is of order 

m th = ^ml + 0{XT 2 ) ~ VXT (2.41) 

when T > m / \/X. 

The imaginary part of the self-energy receives an 0(g 2 T 2 ) contribution from the one- 
loop diagram 0c, but this contribution vanishes for on-shell external momenta since an 



7 This estimate is for the external momentum of 0{T). For a soft external momentum, the diagram ||c is 
0(A 3 / 2 T 2 ). However, this is still 0{y/\) smaller than the diagrams Fig. ||a, b. 
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on-shell excitation of mass m t h cannot decay into two on-shell excitations with the same 
mass. (Nor can an on-shell excitation absorb the momentum of a thermal excitation and 
remain on-shell.) Hence, the dominant contribution to the on-shell imaginary part of the 
self-energy comes from the two-loop diagrams shown in Fig. |3|. At high temperature, these 
two-loop diagrams produce a 0(X 2 T 2 ) imaginary part of the self-energy.^ Consequently, 
at high temperature, the thermal width Fk, as defined in Eq. (|2.38 ), is 0(X 2 T) for hard 
(compared to m t h) external on-shell momenta, and 0(X 3 ^ 2 T) for soft on-shell momenta. 



III. CLASSIFICATION OF DIAGRAMS 

A. Near on-shell singularities of cut diagrams 

Diagrams contributing to the spectral density of the stress tensor correlations function 
have two external vertices each of which connect to at least two propagators. For example, 
the shear viscosity requires evaluating the correlation function 

OW = (TTlmTT 1 " 1 ) , (3.1) 

where the traceless stress tensor, 

TTi m = di<fi d m <p - \5i m d k (j) d k <p , (3.2) 

is quadratic in the scalar field. Naively, one would expect the dominant contribution to 
come from the one-loop diagram shown in Fig. f|. However, in the zero momentum, small 
frequency limit, a finite temperature cut diagram such as this one can contain pairs of lines 
sharing the same loop momenta. As explained below, a near on-shell singularity appears 



Diagrams 3b and 3c for soft external on-shell momenta are 0{g 4 T 2 /mf h ), which is smaller than C(A 2 T 2 ) 
by a factor of 0(m\jrn\^)- Diagram 3d in the same limit is 0(Xg 2 T 2 /m 2 h ), which is smaller than 0(A 2 T 2 ) 
by a factor of 0(rn^/m 2 ^. For hard external on-shell momenta of 0(T), diagram 3a is strictly 0(A 2 T 2 ), 
while diagram 3b is 0(g 4 T/mj? h )(< 0(A 2 T 2 )) due to near-collinear divergences cut-off by the mass, and 
diagrams 3c and 3d are 0{g i /T 2 ) and 0(A<7 2 ), respectively. The explicit evaluation of diagram ||a, at zero 
external momentum, can be found in appendix 0. 
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FIG. 4. A typical one-loop cut diagram for the calculation of a Wightman function. The black dot at 
each end represents the external bilinear operator. 

wherever there is a product of two equal-momentum propagators. Since the thermal width 
that regulates these on-shell singularities is 0(\ 2 ), the size of a diagram is no longer given 
simply by the number of explicit interaction vertices.^ 

The infrared behavior of a cut diagram at non-zero temperature is more singular than 
at zero temperature. At zero temperature, lines in a diagram sharing the same loop 
momentum do not cause on-shell singularities because the poles in the frequency plane all 
reside on one side of the contour. However, at non-zero temperature, a propagator has 
poles on both sides of the contour, as can be seen in Eq. (|2.40aj) . Hence, products of free 
propagators sharing the same loop momentum contain poles pinching the contour, and 
thus produce an on-shell singularity. 

Inclusion of the finite thermal width, as in Eq. ( |2.33| ) and Eq. ( |2.40a| ), regulates these on- 
shell singularities. The effect of these cut-off singularities may be illustrated by analyzing 
the would-be divergent part of the product of two propagators G{k) G(k+5). This product 
represents, for example, the two lines connected to the external vertex on the right side in 



9 In addition to the near on-shell singularities regulated by the thermal width, the soft and collinear 
singularities regulated by the thermal mass must be also considered at high temperatures. Fortunately, 
these soft and collinear singularities turn out not to affect the power counting in A presented in this section. 
Consequently, discussion of this point is deferred until appendix ||]. 
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FIG. 5. Nearly pinching poles in the product G(k) 2 . The heavy line along the real axis represents the 
integration contour. 



Fig. |l] if the small external momentum leaving the vertex is 5. 

As explained earlier, the propagator G(k) has poles at k° = ±(E k —iT k ), and k° = 
±(E k +iT k ). Hence, when 8 — > 0, the product G(k) G(k+8) contains poles separated by 
±iT k on opposite sides of the frequency contour, as illustrated in Fig. [|. When the fre- 
quency integration is carried out, the contribution from these nearly pinching poles is 
0(1/T k ) = 0(1/X 2 ). Exactly the same argument applies to the case of two cut prop- 
agators, or the product of cut and uncut propagators. Hence, the product of any two 
equal-momentum propagators will contain nearly pinching poles. Consequently, a diagram 
with m explicit interaction vertices and n pairs of equal-momentum lines is potentially 
0(X m /T n ) = 0(\ m ~ 2n ). The naive expectation of one-loop dominance is not justified 
when (2n—m) > 2. 

The physical origin of the near infrared divergences caused by nearly pinching poles at 
non-zero temperature can be traced to the existence of on-shell thermal excitations. When 
a small momentum is introduced by an external operator, an on-shell thermal excitation 
can absorb the external momentum and become slightly off-shell. The slightly off-shell 
particle may propagate a long time before it discharges the excess momentum and returns 
to the thermal distribution. Indefinite propagation of a stable on-shell excitation causes a 
divergence, since the amplitude is proportional to the infinite propagation time |TJ]]. But at 
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finite temperature, excitations cannot propagate indefinitely through the thermal medium 
without suffering collisions with other excitations. Hence, there are no stable excitations 
at non-zero temperature. If an excitation with momentum k undergoes collisions at an 
average rate 1/Vfc, the contribution of that mode will be proportional to rfc, or the inverse 
of the width IV 

This may easily be seen explicitly in the product G(k) G(k+5) which contains the 
(nearly) singular piece, 

x(l + 0(A 2 )). (3.3) 

Here the subscript "pp" indicates the pinching pole contribution. The spectral density 
with a Bose factor n(Ek) in Eq. (|3.3|) may be interpreted as available phase space of the 
initial thermal particle. The rest may be interpreted as the Bose-enhanced amplitude for 
propagation of a particle after it has absorbed the soft momentum. When the thermal 
width is small compared to the average thermal energy, the single particle spectral density 
(c.f. Eq. ( |2.37| )) becomes sharply peaked near k° = ±£V Near these peaks, the denomi- 
nator in Eq. ( |3.3j ) becomes O(EkTk)- Hence, the contribution of (G(k) G(k+5)^j contains 
an 0(1/Tk) factor. 

B. Classification 

To simplify the presentation, the classification of A0 4 diagrams will be examined first. 
The effect of adding an additional g<p 3 interaction will be discussed afterwards. 

The classification of the diagrams is fairly straightforward. One only has to count the 
number of explicit interaction vertices in the diagram plus the number of equal-momentum 
pairs of lines as the external 4-momentum goes to zero. Since the thermal lifetime in A0 4 
theory is 0(1/ A 2 ), a finite temperature cut diagram with m interaction vertices and n two- 
particle intermediate states contributes at 0(X m ~ 2n ). For example, the one-loop diagram 
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FIG. 6. The first few chain diagrams in A^> 4 theory. Again, the black circles at each end represent 
bilinear external operators. 

in Fig. [| has a single pair of lines with coincident momenta in zero external momentum 
and frequency limit. When cut, one line effectively forces the other line on-shell, and the 
contribution of the one-loop diagram in the zero momentum limit is 0(1/ A 2 ). 

To determine what diagrams dominate in the calculation of a bilinear operator spectral 
density, one must examine which processes can scatter two particle intermediate states 
into two particle intermediate states. The minimal way of producing a two particle state 
from another two particle state is via a single elementary scattering. Diagrams in A0 4 
theory that consist of only these processes will be called "chain" diagrams. As illustrated 
in Fig. ^, a chain diagram consists of a series of one-loop bubbles. 

Adding each bubble to the chain introduces one additional factor of A from the interac- 
tion vertex and two inverse powers of A from the (nearly) pinching poles of the new bubble. 
Since the lowest order (one-loop) diagram is (9(1/A 2 ), a chain diagram with n bubbles is 
potentially (9(l/A 1+n ). This suggests that the most significant contribution with a given 
number of interaction vertices would come from such chain diagrams. However, the contri- 
bution of each added bubble actually lacks a pinching pole contribution. Consequently, as 
will be shown shortly in section |IIID| , the net contribution of chain diagrams is to modify 
the contribution of the external vertex by a term of 0(XT 2 ). For the bulk viscosity, this 
correction is not negligible since an insertion of V — V—v 2 s (where v 2 is the speed of sound) 
produces an 0(XT 2 ) factor for typical loop momenta of 0(T), as shown in section |IVE| . 



For the shear viscosity, chain diagrams do not contribute at all since the angular integration 
over a single insertion of iri m , (kik m —^5i m k 2 ), vanishes due to rotational invariance. 
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k p 

FIG. 7. A diagrammatic representation of momentum exchange between two lines via one-loop process 
in the A</> theory. When all momenta are on-shell, this process can be interpreted as a second order 
scattering that causes a transition between different two particle states. 

The next most efficient way of causing a transition between different two particle states 
in A0 4 theory is via a second order elastic scattering involving a spectator particle in the 
thermal medium, as illustrated in Fig. In this case, momentum is exchanged between two 
lines via a one-loop process as shown in the first diagram in Fig. |7]. When all momenta are 
on-shell, this process may be interpreted as a second order scattering involving a physical 
thermal particle with momentum I that causes a transition between two particle state 
with a common momentum k and two particle state with a common momentum p. A 
diagrammatic representation of this process is shown in the second diagram in Fig. 0. 

Diagrams in A0 4 theory consisting entirely of two parallel lines exchanging momenta via 
such one-loop diagrams will be called "ladder" diagrams, and are are illustrated in Fig. |8]. 
The one-loop sub-diagrams connecting the other two lines are the "rungs" of the ladder. 
All ladder diagrams contribute at the same order as the one-loop diagram (i.e., Oil/ A 2 )) 
since each rung adds two more factors of A and one more 0(1/ X 2 ) lifetime. Therefore, 
all ladder diagrams must be summed to evaluate the transport coefficients correctly. The 
explicit forms of these ladder diagrams will be examined more closely when the summation 



of all ladder diagrams is discussed in section [TV 



The presence of an additional cubic interaction generates one additional "chain" dia- 
gram and a set of simple "ladder" diagrams whose contribution potentially grows as more 
loops are added. The only "chain" diagram with only cubic interactions is the two loop 
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FIG. 8. The planar ladder diagram with N rungs in \<p 4 theory. The black dot at each end represents 
an insertion of an external operator. 



FIG. 9. The two loop chain diagram in a scalar gift 3 theory. 

diagram illustrated in Fig |9|. Other possible "chain" diagrams with more than two bub- 
bles connected by single lines do not appear because they are a part of the resummed 
propagator. Again, for the shear viscosity, the two-loop diagram vanishes due to rota- 
tional invariance. For the bulk viscosity, as shown in section D fD|, the contribution of this 



two-loop diagram is also to modify contribution of the V vertex by a term of 0(XT 2 ) in 
addition to the modification from summing up A0 4 chain diagrams. The set of diagrams 
that may potentially grow with the increasing number of loops is the set of g<p 3 "ladder" 
diagrams with straight rungs, shown in Fig. |10[ Recall that g = 0(yXm p i l y B ). Hence, 



superficially a ladder diagram with n straight rungs could be 0{1/X n+2 ) since there are 
ri+l factors of 1/A 2 coming from the n+1 pairs of equal momentum lines, and 2n factors 
of g (or equivalently, n factors of A) from the explicit interaction vertices. However, each 
straight rung actually contributes an 0(g 4 ) suppression rather than 0(g 2 ) suppression, 
and hence all ladder diagrams with straight rungs can contribute at C(l/A 2 ), the same as 
the one-loop diagram. 

To understand this suppression, first consider a ladder diagram with the cut running 
through all the straight rungs. When all loop momenta flowing through the side rails are 
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FIG. 10. A ladder diagram with N straight rungs in a scalar g<f> 3 theory. 



forced on-shell by the pinching poles, the momenta flowing through the straight rungs are 
necessarily highly off-shell. Each cut rung contributes a factor of the spectral density 

Ui+k) 2 + m 2 th \ 2 

where { 1; £ 2 are the on-shell 4-momenta flowing through the side rails sandwiching the rung. 
Recall that the imaginary part of the self-energy at an off-shell momentum is 0(g 2 ) = 
0(\m 2 hys ). Hence, when the denominator is 0(m A hys ), a cut rung is 0(g A /Wp hys ), or 
0{\ 2 ). At temperatures comparable or smaller than the physical mass, (T < wiphys), the 
denominator in Eq. ( |3.4j ) is 0(m A hys ) since the typical size of loop momenta is 0(T). 
Consequently, all ladder diagrams with straight cut rungs can contribute at 0{l/\ 2 ) when 
T < friphys- At T = C(m p h ys /"\/A), the denominator in Eq. ( |3.4| ) can be C(mp hys ) = 
0(mf h ) when the small loop momentum contribution cannot be ignored, which is the case 
when calculating the bulk viscosity. At much higher temperatures (T > m phys /A), the 
contribution of a cut rung is at most 0(g 4 / m th) = 0{X A ). Hence, the contribution of a 
ladder diagram containing such rungs may be ignored compared to the contribution of the 
one-loop diagram. 

This additional suppression would appear to be absent when there are uncut rungs. 
This is correct for individual diagrams with uncut rungs. However, as shown in the next 
section, the real part of a rung cancels in the pinching pole approximation when all the 
cut diagrams associated with one original Feynman diagram are summed. Hence, after 
summation over all possible cuts, any straight rung may be regarded as 0(X 2 ). 
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FIG. 11. One-, two- and three-loop diagrams contributing at 0(l/g 4 ). The thick lined sub-diagram 
are all 0(g 4 ) "rungs". Cut lines are not explicitly drawn. 

The key result for the above estimate is that when the loop frequency integrations 
are carried out, the contribution of the sub-diagram sandwiched between pinching pole 
side rails (in this case, the straight rung) can be 0(X 2 ) for T < m p hys/vA. Note that the 
sandwiched sub-diagram need not be restricted to the straight rung for the above estimate 



to hold. Substituting a straight rung with any of the other "rungs" shown in Fig. [TT 
would work just as well, since they all can be 0(X 2 ) when T < m p h ys /"\/A without further 
suppression. 

One important complication is that, for straight g<ft 3 ladders, it is not sufficient to re- 
place the product of propagators representing the side rails by their pinching pole part. 
The non-pinching pole part, G{k)G{k+5) — (G(k)G(k+S)) pp , can also generate leading or- 



der contributions. Specifically, consider the box diagram shown in Fig. [12]. When the 
frequency integration is carried out, the residue of the pinching poles contained in the side 
rail propagators is 0(X 2 ) due to four explicit factors of g from the interaction vertices, 
and one 0(1/ X 2 ) thermal lifetime compensated by two O(X) cut propagators at off-shell 
momenta. This is not the only 0(X 2 ) contribution contained in the box diagram. Putting 
the two cut rungs on-shell also produces an 0(g A ) = 0(X 2 ) contribution since no near- 
divergence cut propagator modifies the explicit factor of g A . It will be convenient to regard 
the 0(X 2 ) non-pinching pole contribution of the box diagram as another elementary "rung" 
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FIG. 12. A box diagram which is a part of a ladder diagram with straight rungs. The heavy lines 
implies that the corresponding momenta are put on-shcll. 



which may be sandwiched between two pinching pole side rails. A more detailed examina- 
tion of the non-pinching pole contribution from the box diagram is contained in section [TV| 
where the summation of all ladder diagrams is discussed. 



C. Higher order rungs in the calculation of the bulk viscosity 

There are other higher-order "rungs" corresponding to processes more complicated than 



those shown in Fig. [IT]. The processes corresponding to these "rungs" contain more ele- 
mentary scatterings than the rungs in Fig. [TTJ without the compensating pinching poles, 
and are sub-dominant as long as individual diagrams are compared. However, when an infi- 
nite number of diagrams are summed, the next order diagrams cannot be simply discarded 
without further analysis of the convergence of the sum of the leading order diagrams. 

For the shear viscosity calculation, no convergence problem arises. However, for the 
bulk viscosity calculation, the sum of the leading order part of the ladder diagrams diverges 
as shortly shown in section [TV]. However, this is not a failure of the theory. As explained 
in section [Tj~B| , the bulk viscosity calculation must involve number-changing inelastic scat- 
tering processes. The leading order part of the simple ladder diagrams contains only the 
elastic scattering processes. Hence, it is no surprise that they cannot produce the correct 
leading order bulk viscosity. 

To calculate the leading order bulk viscosity, the next-to-leading order diagrams 
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FIG. 13. Typical g(jf+\<j) theory 0{\ g ) rungs containing 2-3 scattering processes. 



containing number- changing scattering processes must be included. The lowest order num- 
ber changing process (hence the shortest relaxation time) in the g0 3 +A0 4 theory is 0(\g) 
2-3 scatterings. A few of such "rungs" containing these processes are illustrated in Fig. [13]. 
Other 0(X 2 g 2 ) rungs can be obtained by attaching one more line to the rungs in Fig. [11 
in all possible ways consistent with the theory. Diagrams containing these rungs must be 
included in the bulk viscosity calculation in the g<j) 3 +X(j) 4 theory. 

For the pure A0 4 theory, the lowest order number changing process is 0{\ 2 ). The 0(X 4 ) 
rungs corresponding to these processes can be obtained by attaching two more lines to the 
rungs in Fig. [II] in all possible ways consistent with the \<j) A theory. 

The rest of this section completes the classification of diagrams by showing how the 
chain diagrams modify the external vertex contribution. 



D. Chain diagrams 

Once again, for the sake of simplicity, A</> 4 diagrams are examined first. The analysis 
of the two-loop g4> 3 chain diagram, diagrams with mixed A0 4 and g<j) 3 bubbles, and the 
examination of chain diagrams with more complicated bubbles will follow. For a given 
number of interaction vertices, chain diagrams in A0 4 theory, such as those in Fig. |6|, 
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FIG. 14. Cut two- loop chain diagrams in the \(f) 4 theory contributing to the Wightman function a AA - 

contains the greatest number of pairs of the lines sharing the same loop momentum. A 
chain diagram with n bubbles is potentially 0(1/X n+1 ) because there are n 0(1/ A 2 ) thermal 
lifetimes and (n—1) explicit factors of A from the interaction vertices. However, this is a 
severe over-estimate since the actual contribution of an added bubble lacks a pinching 
pole contribution. This is because (a) the discontinuity of a bubble vanishes in the zero 
external 4-momentum limit, and (b) the real part of a bubble does not contain pinching 
poles. For example, consider the two- loop chain diagrams, depicted in Fig. contributing 
to the calculation of the Wightman function cr AA (q). Here, the external operator A may 
be any component of the stress-energy tensor, and is assumed to be even under a CPT 
transformation. 

The cut bubble is given by 

iH 

(2~y- 

and the uncut bubble in the unshaded region is 

dH 



L A {q) = \ J tSi J a(U-0 5(0 S(q-l) , 



(3.5) 



C A (q) 



T 



(27T)' 



I A (l,q-l)G(l)G(q-l) 



(3.6) 



where q is the external 4-momentum, and Ia{1, denotes the (polynomial) contribution 
from the external operator in such a manner that the contribution of <f) 2 1^(1, q—l) = 1. 

Since the operator A is even under a CPT transformation, at zero external momentum 
Ia{1, — I) is a real, even function of the loop momentum I. Hence, the sum of the two-loop 
chain diagrams in the zero external 4-momentum limit is 

hm lim LSS oop) (g) = 4XL A (0) (C A (0) + C A (0)*) 

q"— >o q^u 
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= 8\L A (0)ReC A (0) . (3.7) 

In the same limit, the cut one-loop bubble La(0) is 0(1/ A 2 ) as before. To see that 
AReC*A(0) does not exceed 0(1), consider the following explicit form of the real part 
of an uncut bubble at zero momentum, 

ReC A (0) = i| -£Ll A (l,-t) Re (-iG\l)) 

=y^m-i)cotHiw ( [l2+m i +m]2 - [P+ ^ +m * ? ) ■ 

(3.8) 

Here, /^(i, — ~ Z 2 if A is a stress-energy tensor. Since the integrand does not contain 
pinching poles (i.e., there is no l/|Z 2 +m 2 h -|-E(/)| 2 term) no large lifetime factor appears 
when the frequency integration is performed. In appendix [F|, the real part of the uncut 
one-loop diagram is shown to be 

ReC A (0) = 0(T 2 ) , (3.9) 

using the fact that the integrand is appreciable only when I is nearly on-shell. 

Individual higher order chain diagrams with more one-loop bubbles strung together 
may be analyzed in a similar manner. However, since chain diagrams form a geometric 
series, it is also straightforward to sum all cut chain diagrams with one-loop bubbles and 
examine the result of the summation. Of course, one can also perform the geometric sum 
first in imaginary-time, and then take the discontinuity of the result of summation. 

The summation of cut chain diagrams with one-loop bubbles is fairly simple. The only 
subtleties come from the cuts involved and the fact that there is an external operator at 
each end of a cut diagram. Due to the cuts, the equation for the resummed chain is a 
matrix equation instead of a single component linear equation. The presence of external 
operators implies that the bubbles at each end are not equivalent to the other bubbles. 

Since no additional difficulties than those already present in the two-loop calculation 
appear, performing the actual summation of the cut chain diagram is deferred to ap- 
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pendix |F|. The result of the summation of all chain diagrams with one-loop bubbles is 
shown in appendix [F] to be 

Km rr^(n) - AT (0) + 8AReC ^°) j (n) + 4(AReC A (0)) 2 
hmo^ (q) - 4 WO) + 1 _ XRjeCo{0) L ^) + (1 - AReC (0)) 2 

= 4^(0), (3.10) 
where the finite temperature optical theorem 

Im C (q) = ~ (L (q) + L (-q)) , (3.11) 



is used to simplify the result. (The optical theorem (|3.11|) can be easily proven from 



Eq. (|3.5|) and Eq. ( |3.6|) .) Here, cr A h \ an (q) denotes the contribution of these chain diagrams 
to the correlation function <taa(<?), Laa(q) corresponds to the contribution of the one- loop 
diagram with the external operator A at both ends, and L A (q) denotes the contribution 
of the one- loop diagram with A at one end. L (q) and Co(q) are the cut and the uncut 
bubbles with I a = 1. The modified one- loop contribution Laa contains the (modified) 
vertex contribution 

= I A (l, -I) + XReC A {0) x (l + 0(V\)) , (3.12) 

where the estimate AReC(O) = 0(V\) is used. This estimate of AReC(O) is justified in 
appendix 0. For the operator V = V—v^e required for the bulk viscosity, Ip = 0(XT 2 ) for a 
typical 0(T) loop momentum, as shown in section [IV E . In the same section, the additional 



term AReC-p(O) is also shown to be 0(XT 2 ). Hence, the correction term AReC^O) in 
Eq. (|3.12|) cannot be simply ignored. For the shear viscosity, Re C n (0) vanishes due to 
rotational invariance. Hence, no modification is needed in that case. 

When cubic interactions are added, the "chain" diagrams also include the two-loop 
diagram shown in Fig. ^| where each bubble in the diagram now may be regarded as the 
sum of all A0 4 chain diagrams. The sum of all chain diagrams in the <?0 3 +A(/> 4 theory is 
given by the sum of the 4 chain result L AA ( |3.10|) and this two-loop diagram. As shown 
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in appendix [F], a straightforward application of the cutting rules yields the sum of all chain 
diagrams as 

lim a^ chain (g) = 4L^(0)-8^-ReC i (0)L i (0)+4^ (Re ^(O)) 2 L (0) 

= AL AA (0), (3.13) 

where L AA (0) contains the modified vertex contribution I A given by 

I A (l,-l) = I A (l,-l) - ^ReC A (0) 

= I A (l, -I) + (\-JL.) ReC A (0) x (l + 0(V\)) . (3.14) 

m th 

More complicated chain diagrams can be produced by including more complicated 
bubbles such as ladder diagrams. For these more complicated "bubbles" , exactly the same 
argument given above will also apply provided that the generalized finite temperature 
optical theorem (c.f. Eq. (|3.11| )) 

ImC* bubble (g) = -- (L bubb i e (g) + L bubblc (-g)) (3.15) 

holds for each bubble. However, unlike the imaginary part, the real part of higher order 
contributions, including ladder diagrams, to the bubble C A are suppressed compared to the 
real part of the one loop contribution. Hence they can be safely ignored. The generalized 
optical theorem fl3.15|) can be inferred from the works of Kobes and Semenoff [15) and will 
not be further discussed in this paper. 



IV. SUMMATION OF LADDER DIAGRAMS 

As explained in section [11 Bj , calculations of the shear and the bulk viscosities require 
different set of diagrams. More specifically, to evaluate the leading order shear viscosity, 
summation of only the leading order ladder diagrams is needed. Whereas, to evaluate the 
bulk viscosity, as shown in this section, C(A 4 ) rungs must also be included. In this section, 
the leading order ladder summation for the shear viscosity is examined first. More 
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FIG. 15. The diagrammatic representation of the integral equation for the effective vertex function 
T>T r {k,q—k). This corresponds to Eq. (4.1). 



complicated analysis of summing the higher order contributions for the bulk viscosity 
follows. The results presented in this section are valid for all temperatures. 

A. Ladder summation for the shear viscosity calculation in A</> 4 theory 

Cut ladder diagrams form a geometric series, and can be resummed by introducing 
a suitable effective vertex. Due to the various possible routings of the cut, the integral 
equation will involve a 4x4 matrix valued kernel. Hence, it is convenient to introduce an 
effective vertex T> n (k,q—k) which is a 4-component column vector. The subscript tt is a 
label for a component of the traceless part of the stress tensor. The resummed effective 
vertex satisfies the following linear integral equation 

V w (k,q-k) = T„{k,q-k) + J -^M{k- P )F{ P ,q-p)V^q-p) (4.1) 

illustrated in Fig. ITR Here, I n (k,q—k) is an inhomogeneous term representing the action 



of the bilinear operator ft including the contribution of chain diagrams, Ai(k-p) is a 

33 



4x4 matrix representing the rungs of the ladder which consist of cut and uncut one- 
loop diagrams, and !F(p, q—p) is a matrix representing the side rails of the ladder and 



consists of products of propagators. As shown in Fig. [15], the first component of T> n (k, q—k) 
corresponds to the effective vertex where momenta k and q—k enter vertices in the unshaded 
region. For the second component, k and q—k enters vertices in the shaded region. In 
the third component, the momentum k enters a vertex in the unshaded region while the 
momentum q—k enters a vertex in the shaded region. The last component of T> n differs 
from the third component by changing k to q—k, and vice versa. 

In a more symbolic form, the above equation can be compactly rewritten as 

|JA = (1-/C)|2\>, (4.2) 

with the identification of the "ladder kernel" 

K = MT. (4.3) 

As is evident in Fig. [T5|, only the first component of the inhomogeneous term I n is non-zero 
and given by ikik m — 5; m k 2 /3). Explicit expressions for M. and T are given in appendix [B]. 
Note that all quantities depend on the external 4-momentum q = (u, q). 

The integral equation \T n ) = (1— K.)\T> W ) will be solvable only if any (left) zero modes of 
the kernel (1— K) are orthogonal to the inhomogeneous term \I n ). The operator (1— JC) does 
have four zero modes in the zero momentum, zero frequency limit. These four zero modes, 
denoted (V M |, are related to insertions of energy-momentum density T M °, and the existence 
of these zero modes is a direct consequence of energy-momentum conservation. Explicit 
forms of these zero modes are shown in appendix y. Reassuringly, is orthogonal to 
the zero modes; this is also verified in appendix |C|. 

In terms of the resummed vertex \T> n ), the Wightman function of a pair of 7T; m is simply 

a n7T (q) = 2 (^l-FlAr) x (1 + 0{\ 2 )) , (4.4) 

where z % represents the action of the operator 7T; m in the same way X w represents the 
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action of the operator 7T; m .|^| The only difference between and X n is that z v has its only 
non-zero component in the second slot while T v is non-zero in the first slot. The overall 
normalization constant 2 is chosen for convenience. 
In this notation, the shear viscosity 77 is given by 

rj = ^ lim lim .(z«\F\V«) x (1 + 0(A)) . (4.5) 

10 w-»0 |q|->0 

From now on, the external 4-momentum may simply be set to zero. 

In the limit of vanishing external momentum, the leading weak coupling behavior is 
generated by the (nearly) pinching pole contribution to the p° frequency integral. Hence, 
portions of the side rail matrix T which do not contain pinching poles may be neglected. 
Examination of appendix [B] together with the explicit form of the cut and uncut single 
particle propagators shows that the leading order part of the remaining pinching pole part 

isQ 

^ PP (P, -V) = w(p)u T (p) x [l+n(E p )\ n(E p ) sgn(p°) 27r5(p 2 -E 2 p ) / £/(p) (4.6) 

where 

w T {p) ee (1, 1, (l+e"^)/2, (l+e^)/2) , (4.7) 
u T {p) ee (1, 1, (l+e^)/2, (l+e- p °^)/2) . (4.8) 

The leading order kernel is given by 



10 The appropriate inner product is defined by 

(f\g) = J ^/(p)t s(p ). 

11 Actually, there is one more part of T that contains pinching poles in the zero momentum, zero frequency 
limit. However, this part, denoted hj T in appendix^, does not contribute to the leading order calculation 
for the following reason. First, hj T is orthogonal to the vertex parts since z T A h — 0, and j T lB — 0. Second, 
hj T is orthogonal to the rest of T in the sense that if f = T—hj T then j T M.T = and TM.h = 0. Hence, 
the hj T part does not affect the contributions of the ladder diagrams in a^l = z\ {? 'M) n TIb ■ 
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/C pp EE Mo^pp , (4.9) 

where M.q contains one- loop rungs evaluated with free propagators, and the self-energy 
S/ in jF p p contains only the contribution of the two-loop diagram calculated with the free 
propagators. Since the factors of coupling constants from Ai and £/ cancel each other, 
/Cpp is independent of A except for those contained in the thermal mass. Note that dropping 
non-pinching pole contributions reduces T to a rank one matrix. This allows one to greatly 
simplify the equation. 

For the change in the solution of the integral equation (O) caused by the replacement 



of /C by /Cpp to be sub-leading in A, the inhomogeneous term must be orthogonal to the 
(left) zero modes of (1— /C pp ) as well as orthogonal to the original zero modes of (1— /C). 
Otherwise the reduced integral equation \X K ) = (l~JC pp )\V n ) would be singular implying 
that the neglected part of /C could not be negligible. The issue of zero modes of (1— /C pp ) 
does not arise when considering the size of an individual diagram as in section |T|, but 
rather reflects the convergence (or lack thereof) of the infinite series of ladder diagrams. 
Suppose the inhomogeneous term X n had a non-zero projection onto a zero mode y. Then 
p P ) n y — V, and all ladder diagrams would contain an identical 0{\/\ 2 ) piece, z^T pp y, as 
a part of their pinching pole contribution. The infinite number of such terms would make 
the sum diverge. Hence, to produce a finite result, the inhomogeneous term must satisfy 

(VslXr) = , (4.10) 

where (b^ )5 \ denotes the five zero modes of (1— /C pp ) whose explicit forms are 

fo» ee (6» = p, [l+n(p )] S hee (-P) u T (p) , (4.11) 

and 

6 5 (p) = (h\p) = sgn(p°) [l+n(p )] S ivcc (-p) u T (p) . (4.12) 

Here, b^s corresponds to the 4- momentum conservation, and the additional 65 corresponds 
to the particle number conservation. Of course the theory does not preserve the number 
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of particles. However, the number changing scatterings are 0(X A ), and hence, do not 
contribute at the leading order. 

As a simple consequence of rotational invariance, the traceless stress operator involved 
in the calculation of the shear viscosity does satisfy Eq. ( |4.10| ). When t>i(k) is applied to 
I n (k), it vanishes since rotational invariance requires that any rank 3 spatial tensor with 
two symmetric indices be a combination of kik\k m and ki5i m . Applying bo(k) or 65 (k) again 
results in zero because the angular integration over X^{k) vanishes. 



The well-posed integral equation (|4.2|) = (1— /C pp )|2\), can now be reduced, since 
the pinching pole kernel J-" pp ( |4.6|) is a rank one matrix, by applying u T to both sides of the 



vector equation. The resulting linear integral equation is (dropping sub-leading corrections 
suppressed by 0(A)), 

I n (k) = D w (k) -j-0ri K pp (k,p)n(p°) S ivcc (p) f^j , (4-13) 

where I n (k) is the first non-zero entry of (k\X n ), the reduced effective vertex is 

D n (p)=u T (p)(p\V n ) , (4.14) 

and the reduced integral kernel is 

K pp (k,p) =u T (k)M {k-p)w{p) 

= I(l_ e -*°/>)i (fc_p)( e pP/»_l). (4.15) 

Here, the explicit form of the free cut particle propagator, 

S {ree (p) = [l+n(p )} sgn(p°) 2nS(p 2 -E 2 p ) , (4.16) 

is used, and L (k—p) is the cut rung given (in A0 4 theory,) by 

L (k-p) EE y J -j^- S ivcc (l+k-p) SfreeH) . (4.17) 

Note that K pp contains no reference to the real part of the uncut rung. When u 7 A4qw is 
calculated, the real part of the uncut rung cancels. Eq. ( |4.15|) is obtained by expressing 
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the remaining imaginary part of the uncut rung in terms of Lo(k-p) with the help of the 
optical theorem ( |3.11| ). 



Due to the delta function present in the kernel, p is an on-shell momentum. Also, since 
the leading weak coupling behavior of Wightman function is given by 

= 2 / ^z T n (k)w(k)n(k°) SUk) ^ , (4.18) 

the final integral over k will be also restricted to on-shell momenta. Hence, the reduced 
integral equation ( |4.13|) need be solved for only on-shell momenta. 



To summarize, after summing all ladder diagrams in A0 4 theory, the loop frequency 
integrals may be performed and the leading weak coupling behavior extracted from the 
pinching pole contribution. The resulting linear integral equation for the effective vertex 
reduces to a single component equation given explicitly by 

KH) = D Ak) - I ^ L 0ik - P) g±gj| W) , (4.19) 

where k is an on-shell momentum. 

The cut rung L (k—p) is easily shown to satisfy 

L (k-p) = e (fc °- p0)/3 L (p-k) . (4.20) 

Also, £/(p) is an odd function of p°. Consequently, D w (—p) satisfies the same equation as 
does D n (p), provided I n (k) is an even function of k. Hence, if ^(k) is an even function 
of k, so is the solution D n (k). Since the energy-momentum tensor is even under CPT, 
the inhomogeneous terms for both the shear and bulk viscosities are even functions of the 
4-momentum. 

In terms of the solutions of the reduced integral equation ( |4.19| ), the shear viscosity is 
V = i; / ^z:(k)w(k)n(k )SUk) D ' J '' } 



10 J (2tt) 4 n " ' v ' v ' ™ v ; E/(fc) 

iJ^rMWs,,^)^ (4. 21 ; 



10 J (2vr) 4 nv ' v >~™^> E 7 (fc) 
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neglecting sub- leading contributions suppressed by additional powers of A. 

For the future use, we define the inner product of two functions of on-shell momentum 

as 



(f\g) 



dH 



n(i°)s iTee (i)f(iyg(i) ■ 



(27r) 4 E 7 (Z) 

In terms of this definition, the integral equation ( j4.13|) can be expressed as 



(4.22) 



\M = (l-K m )\D 



(4.23) 



whose five zero modes are 



and 



(b,\p)=b ll (p)=p ii Z I (p) 



(h\p) = hip) = sgn(p°) Ej(p) 



(4.24) 



(4.25) 



B. Ladder summation for the shear viscosity calculation with an additional gcf) 3 

interaction 

To start, consider "simple" ladder diagrams only containing the straight single line 



rungs, as illustrated in Fig. [10]. After summing these diagrams, including the contribution 
of the other required rungs will be easy. To sum these simple ladder diagrams, one again 
introduces an effective vertex T> n (k,q—k). Before performing any frequency integration, 
the effective vertex satisfies 

V v (k, q-k) = l^k, q-k) + J M linc (k~p) F(p, q-p)V„{p, q-p) , (4.26) 

where the elements of the matrix M.n ne {k—p) are simply cut and uncut single particle 
propagators. Before proceeding with the general analysis, it may be helpful to consider a 
typical example, such as the three-loop diagram in Fig. [Tj| Applying the cutting rules, the 
contribution of this three-loop diagram (with zero external 4-momentum) is 



39 



k—l 




777777 k 



l-k 



FIG. 16. A cut three loop ladder diagram in a scalar g(j> 3 theory. The black dot at each end represents 
an insertion of a bilinear external operator. 



(0) = 2 J — jJL I n (k) S(k) G(k)*L Ml _ hox (P, k) S(-p) G(p)h(p) , (4.27) 



The size of the contributions from the (nearly) pinching poles in the complex k°, p° 
and 1° planes can be estimated as follows. Each pinch generates a thermal lifetime of 
0(1/ X 2 ). There are four 0(\j r Xm p \ iys ) factors of g. When all three loop momenta k, p 
and I are on-shell, the momenta flowing through the two cut propagators are well off-shell. 
An off-shell cut propagator is 0(g 2 ) since it is proportional to the imaginary part of the 
0(g 2 ) = 0(Xm 2 hys ) one-loop self-energy (c.f. Eq. ( |2.33|) ). Hence, when T < m phys , the 
pinching pole contribution is g A x 0(1/ A 6 ) x 0(g A ) = 0(1/ A 2 ), or the same as the lowest- 
order one-loop diagram. Since the leading order shear viscosity is insensitive to small 
momentum contributions, at temperatures much greater than m p h ys , the cubic interaction 
becomes irrelevant compared to the quartic interaction.^ 

Equivalently, when T < m p h ys the pinching pole contribution of the box sub-diagram 
is 0(g s /X 2 ) = 0(g 4 ), not 0(g 4 /X 2 ) = 0(1) as one might have expected if the additional 



12 In contrast, the bulk viscosity is sensitive to small momentum contributions, and therefore the cubic 
interaction becomes negligible at much higher temperature T 3> rn p hy S /i7\. 



where 




(4.28) 



is the one-loop box sub-diagram illustrated in Fig. [L7|. 
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l—k 

FIG. 17. A box diagram which is a part of a ladder diagram with straight rungs. 

suppression from the off-shell self-energy were ignored. Consequently, the non-pinching 
pole contribution of the box diagram, which is also 0{g i ) from the four explicit interaction 
vertices, is equally important as the pinching pole contribution. This complicates the 
treatment of these diagrams. 

The key observation of the above argument is that an off-shell straight cut rung is 
0(g 4 /m,p h y S ) = 0(X 2 ) when T < m p h ys . Hence, the leading weak coupling behavior of the 
non-pinching pole contribution is produced when the two cut propagators in the box are 
both on-shell. Otherwise, there will be additional suppression from the cut propagators 
which will make the contribution smaller than 0(X 2 ). Consequently, the contribution of 
the three-loop ladder diagram in Fig. WQ can be rewritten as 



a (Fi,16) (0) = 

74 7 7 4 

+ 2 / my* mh h{k) Lbox(p ' k)J7 ™ h{p) ' (429) 

where is the (ij) element of the pinching pole side rail matrix jF pp , and the "rungs" 
between the pinching pole side rails are 

£box(p,£) =9 A f G(k—l) G{l-k)* S frcc (-/) S hee (l-k+p) , (4.30) 
representing the 0(X 2 ) non-pinching pole contribution from the cut box sub-diagram 
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FIG. 18. Diagrammatic representation of a cut line rung corresponding to Eq. (4.31). The heavy line 
in the unshaded region represents a retarded propagator Gr, and the heavy line in the shaded region 
represents an advanced propagator Ga- 



(where the resummed cut propagators may be replaced by free cut propagators), and 
L\ine{L-k) = 9 2 



[l+n(f-fe°)] 2 Sj(l-fe) 

\a-kr+mi + m~k)\ 2 



= g 2 [l+n(l°-k )] 2 £/(/-£) G R (l_-k) G A (l_~k) 

= ^ L{l_-k)G R (l_-k) G A {l-k) , (4.31) 

representing a single cut straight rung as shown in Fig. [Tj|. The subscripts R and A here 
indicate the retarded and the advanced propagators (c.f. Eq. (|2.35|) ). 

Applying the same argument above to other routings of the cut in the three-loop ladder, 
it is straightforward to see that the sum of all cut three-loop ladder diagrams has the form, 

a (3-loop) (0) = 

d 4 k d A p 
(2tt) 4 (2tt) / 

where a 4x4 matrix Mh ox (k,p) contains non-pinching pole contributions of cut and un- 
cut box sub-diagrams. The previous expression ( [4. 29| ) is included since, as shown in ap- 
pendix [FJ, the (33) component of M.\\ne{k—l) is Lr mo (Z— k), and the (33) component of 
■Mbox(k,p) is Lh ox (p—k). (Recall also that 2 n = (I n , 0,0,0) and z n = (0,7^,0,0).) Once 
again, merely replacing the side rail matrix T by its pinching pole part jF pp , without 
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+ 2 / ^^^(k)T pp (k)M hox (k,p)^ pp (p)Up) , (4.32) 



changing the rung matrix .Mime, is not sufficient to calculate the the leading weak coupling 
behavior of a simple ladder diagram; one must also include the box sub-diagram rung 

-Mbox- 

To sum all simple ladder diagrams, note that the first term in Eq. (|4.32|) can be regarded 
as the second iteration of the single-line kernel -Mi me .Fp p , and the second term can be 
interpreted as the first iteration of the box kernel A^box^pp- In exactly the same way, it 
is simple to deduce that the leading weak coupling behavior of a simple ladder diagram 
with n straight rungs contains all possible sequences of n—2m factors of single line rungs, 
-Miine^pp, and m factors of box rungs, A^box^pp, for all m < n/2. Every such sequence 
can be interpreted as arising from the iteration of the integral equation 

l n (k)=V 7T (k)- I (M h Uk-p)+M hox (k,p))F pp ( P )V w (p) . (4.33) 

Because the pinching pole kernel JF pp ( |4.6| ) is a rank one matrix, applying u T to both 
sides of Eq. ( |4.33| ) reduces the equation to 

h{k) = Ar® - / 7^ K simplc . laddci (k,p)n(p°) S ircc (p) ^44 > ( 4 - 34 ) 

where 

^simple-ladder (k,p) = ^ ( 1 — e ~ ^ ^ ) -^simple-ladder (k, p) (e P ° ^ — 1) (4.35) 

with 

f dH dH 

^simple-ladder (&,P) = J ^— - ^-yj (2tt) 4 5(/i -l 2 +p-k) Sfree^l) S hcc (~l 2 ) 

x (^G R {k-p) G A (k-p) + ^G R (k-h) G A (k-h) 
+ g -^G R (k+l 2 )G A {k+l 2 )^ . (4.36) 
In obtaining Eq. ( fl.34j ), the following relation for the box rung is used, 

u T (k)M hox (k,p)w(p) = (1-e-^) (M^(k, P )-e k °^M^(k,p)) (e^-l) I2 , (4.37) 
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together with an analogous relation for .Mime and Li me . Verifying equations ( |4.36|) and 
( |4.37 ) is a straightforward but tedious exercise in the application of the generalized optical 



theorem. In short, to prove Eq. (|4.37| ), one must show that all cut box diagram contribu- 
tions in u T AiboxW other than the (44) and (34) contributions are canceled by the imaginary 
part of the uncut box diagram. The final expression (|4.36|) in terms of the retarded and the 



advanced propagators results from the particular combination of (44) and (34) components 
in Eq. (|4.371 ). A sketch of the proof is given in appendix |B[ 



The above summation of simple ladder diagrams illustrates the general principle: To 
determine whether a diagram contributes to the leading weak coupling behavior, first carry 
out the frequency integrations. Then, if the contribution of each sub-diagram sandwiched 
between two pinching pole side rails is 0(X 2 ), the diagram as a whole will contribute 
to the leading weak coupling behavior. Hence, to identify all diagrams in the g0 3 +A0 4 
theory contributing to the leading order behavior, one must identify all 0{\ 2 ) sub-diagrams 
("rungs") which may be sandwiched between two pinching pole side rails. In the simple 
ladder diagram, two of such "rungs" were identified, the single line "rung" and the box 
sub-diagram "rung", illustrated in Fig. |T§| labeled (b) and (c). 



With cubic and quartic interactions, there are a total of 10 different 0(X 2 ) "rungs" 
when T < m p h ys as shown in Fig. These diagrams exhaust all possible 0(X 2 ), 0(\g 2 ), 
0(g A ) rungs in g</> 3 +A</> 4 theory in those temperature range. The particular cuts shown 
in the figure correspond to the (44) components of the rung matrices. Consequently, the 
dominant set of "ladder" diagrams in g0 3 +A0 4 theory may be described as iterations of 
pinching pole side rails T vv and a combined 0(\ 2 ) rung matrix .M m ii, whose components 



are the sum of contributions of all possible cuts of the underlying diagrams of Fig. [19[ The 
summation of all such ladder diagrams is generated by the integral equation 

l w (k) = V n (k) - J A^MmUhp) T pv (jp)V«{ v ) , (4.38) 

where -M m u(/c,p) consists of all cut and uncut 0{\ 2 ) rungs in the g0 3 +A0 4 theory, i.e., 
those in Fig. [H| Once again, applying u T to both sides of Eq. Q4.38j ) reduces the equation 
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(f) (g) (h) (i) 0) 

FIG. 19. The C(A 2 ) "rungs" in the grfP+Xcj) 4 theory. The cuts shown correspond to (44) components 
of the rung matrix. 



to 



I w (k) = D^k) - J K m (k,p) n(p°) S hcc (p) ^ 



DM 

(P) 



(4.39) 



where the previous relation between the rung matrix -Mf u n and its (44) and (34) components 
continues to hold, 



K M i{k,p) = u T (k)M iu \\{k,p)w{p) 

= \(l-e- k0 P)L lull (k,p) (e^-1) 



with 



P 



(4.40) 



(4.41) 



as a consequence of the generalized optical theorem. The proof of Eq. (|4.40| ) is discussed 
in appendix [B]. 



Noting that all the diagrams in Fig. [19] have two on-shell cut lines, a straightforward 
application of the cutting rules yields the sum of all the cut "rungs" , 
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f dH dH 

Lm(t,p) = J j^y j^y S'free(^) S frcc (-l 2 ) (2tt) 4 S(k-l 2+p-k) 

X (y + J Gr ^~S) G ^ k ~^ + J G x(k-h) G A {k-h) 

+ 9 -^G R {k+l 2 )G A {k+l 2 ) 
+ g 4 Re (G R (k+p)G A (k-k)) + g 4 Re (G R {h+p) G A {h-l 2 ] 
+ 4 Re (G R (h-k)G A (h-h] 



- Xg 2 ReG R (h+p) - Xg 2 ReG R (h-k) - Xg 2 ReG R (h-l 2 )j (4.42) 

5 / Si 7^4 ^rccfa) S frcc (-/ 2 ) (2 7 r) 4 5(/ 1 -/ 2 +p-i) 

x A - <j 2 (G^+p) + G^-A;) + G^-^)) % . (4.43) 



Each term in (|4.42|) arises from specific diagrams of Fig. [19|. For example, consider the 
term g 4 Re (G R (h+p) G A (l\— ku . This term corresponds to the "cross" diagram labeled 
as (d) in Fig. and redrawn in Fig. |2(] with momentum labels. The cutting rules produce 

W&P) = M<?±(k,p) - e^Mf±(k,p) 

X (Sfree(Zi) 5frcc(-/ 2 ) G(g+/l) G^-*)* - S^-h-p) S ttee {l X -k) G{1 2 ) G(h)*) 

x (SUfr) S froc (-Z 2 ) G«(p+Z 1 ) G A {h-k)-e^ SfeeHi-p) S^h-h) G R {1 2 ) G A {h) 



(4.44) 

where to obtain the last expression, the relation between propagators 

G(k) = —iG R {k) + S{-k) = -iG A {k) + S{k) , (4.45) 

and the fact that S(h — k) and S(p + h) can be neglected since l\, l 2 ,P,k are all on shell are 
repeatedly used. Noting that the effective vertex D n (p) is an even function of p, the sign 
of p in the second term of Eq. ( |4.44| ) may be changed inside the integral equation ( |4.3<j| ). 
Then changing the label l\ — > h+p, replacing full cut propagators by the free cut 
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—k p k p 

l\—k h 

FIG. 20. The (44) and (34) components of the "cross" sub-diagram rung. 

propagators (and dropping sub- leading corrections suppressed by 0{\/X)) yields, 

d% d% 



L cmss {k,p) =9 A J j^Ti T^h (27r) 4 5(/i-/ 2 +p-i) 



(2tt) 4 (2tt) a 

x S bee (h) 5 frcc (-/ 2 ) Re (G R (h+p) G A (h-k)) . (4.46) 



All other terms are produced similarly. 

The most important point to notice in Eq. ( f4.43|) is that the various terms combine to 
produce the square of a single factor 

TO^ia, AO = X- g 2 (Gnik+p) + G R {l x -k) + G R {l x -l_ 2 )) , (4.47) 

which obviously resembles a tree level two-body "scattering amplitude" . Strictly speaking, 
at non-zero temperature, one cannot define scattering amplitudes since no truly stable 
single particle excitations exist. However the expression ( [4.47[ ) may be regarded as an 
approximate scattering amplitude characterizing the dynamics of the finite temperature 
excitations on time scales short compared to their lifetime. The only difference between the 
result for g0 3 +A0 4 theory, and that for the A0 4 theory, is the replacement of the constant 
tree level scattering amplitude A in the A0 4 theory by the momentum dependent tree 
amplitude T(i 1; p; lz,k). Note that T contains retarded propagators in place of the usual 
time-ordered propagators. At zero temperature, the scattering amplitude can be expressed 
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both in terms of the time ordered Green function, or the retarded Green function jLf]]. At 
non-zero temperature, it is retarded Green function which gives the correct amplitude. 
The arguments of the propagators in Eq. ( |4.47| ) are all combinations of two on-shell 



momenta. A short exercise in kinematics shows that the 4-momentum squared of the sum 
of two on-shell momenta is always less than —4m 2 h , while the 4-momentum squared of the 
difference of two on-shell momenta is always greater than 0. Hence each propagator in 
( [4.471 ) is bounded by l/m 2 h so that 

{TfanhM = 0(\) + 0(g 2 /m 2 h ) = 0(A) , (4.48) 

since g 2 = 0(\m^ hys ) and (^ph ys / m th) — 1- Consequently, the size of the effective ver- 
tex D^ip) with the g0 3 +A0 4 kernel does not differ (by powers of couplings) from the 
solution with only the A0 4 kernel. Furthermore, formula (|4.21|) for the shear viscosity con- 



tinue to hold in the g4> 3 +\4> 4 theory (with the self-energy now given by the full g0 3 +A0 4 
self-energy). Hence, in both theories, the shear viscosity is (9(1/A 2 ), only the numerical 
prefactors will differ. Note also that at temperatures T 3> m p h y s, the typical size of a prop- 
agator is 0(1/T 2 ) and hence the contribution of the cubic interaction g 2 G^ = 0(g 2 /T 2 ) 
may be ignored compared to A. 

For the comparison with the kinetic theory in the next section, it is helpful to note 
that the imaginary part of the on-shell two- loop self-energy E/(fc) (c.f. Fig. |^), used in 
Eq. ( f4.39|) , can be expressed in terms of the same "rungs" of Fig. [19| by joining two of the 



external lines by a cut propagator, 

Mt) = l(l-e- pk °) J LUk-p)S bee (p) . (4.49) 

For instance, the contribution of the A0 4 self-energy diagram can be obtained by attaching 
a cut line to the A0 4 rung labeled (a) in Fig. [19] and dividing by 6. Since the diagram (a) 
has a symmetry factor of 2, this correctly reproduces the overall factor of 1/12 associated 
with this self-energy diagram (6 from the symmetry factor for this two-loop diagram, and 2 
from the relation between the imaginary part of a diagram and the discontinuity). The g<j) 3 
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theory two-loop self-energy diagram labeled (b) in Fig. |3| is similarly obtained by attaching 
a cut line to the rungs labeled (b) and (c) in Fig. [19] and dividing by 6. Since the diagram 
(b) has a symmetry factor of 2, and the box diagram (c) has a symmetry factor of 1, this 
again correctly produces the overall factor of (l+l/2)/6 =1/4 associated with this self- 
energy diagram. All other two-loop self-energies can be reproduced from the rung diagrams 
in a similar manner. 



C. Ladder summation for the bulk viscosity calculation in Xcj) 4 theory 



As explained in section |IIB| , when the volume of a system changes (or equivalently when 
the density changes), number-changing processes are ultimately responsible for restoring 
equilibrium. Hence, the calculation of the bulk viscosity must include the effect of changing 
particle number. The leading order equation 

\Xp) = {l-M in iiF pp )\V P } , (4.50) 

where A4f u \\ and jF pp are given in Eq. ( |4.38| ), contains only the effect of elastic binary 
scatterings, and therefore is not suitable for the calculation of the bulk viscosity. 

Mathematically, the integral equation (|4.50| ) is not well posed because only one of the 
consistency conditions (b \Tp) = and (b 5 \Xp) = can be satisfied by adjusting the value of 
a single free parameter v 2 in V. Since energy-momentum must be conserved, the condition 
\Ip) = must be enforced by choosing 

■n(E k )[l+n(E k )}I v (k) 



Jo 



v2 = (bo\Xp) _ J (2tt) 



7 1 j^y 3 n(E k )[l+n(E k )]I E (k) 

(4.51) 



(dV/dT) , 



(de/BT) 

Here, Ij>(k), and I £ {k) represent the effect of the pressure and the energy density insertions 
including the contribution from chain diagrams (c.f. Eq. (|3.14j) ), and v s is the speed of 
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sound.0 Note that the condition v 2 = v 2 is a familiar result from the Boltzmann equation. 
In section [IV E] , it is shown that with this choice of v 2 , Xp oc m 2 hys . The (leading order) 
bulk viscosity vanishes if the mass is zero since Xp = in that case |J. 

When (9(A 4 ) rungs are included in the ladder kernel, 65 is no longer a zero mode since 
the kernel now contains 2-4 number-changing processes. Since the condition (6 5 |Ip) = 
need no longer be satisfied, the resulting integral equation is well-posed, and hence the 
leading order bulk viscosity can be evaluated as follows. 

The integral equation for the effective vertex including 0(X 4 ) rungs may be written as 

|Xp) = (l-/C cons -5/C ch .)|%) , (4.52) 

where /C cons . is the number conserving part of the kernel, and 5/C c h. is the number changing 
part of the kernel. In terms of the solution \T>p), the leading order bulk viscosity is given 
by 

C = (3(Ip\F\Vp) . (4.53) 

As will be shortly shown, 65 is no longer a zero mode of the kernel due to the number 
changing <5/C c h. part. However, b 5 is still a zero mode of the number conserving part of the 
kernel so that 

<6 5 |(l-/C cons .- ( 5/C ch .) = -{b 5 \81C ch . . (4.54) 

Since this vanishes as A —>■ 0, the kernel (1— /C cons — 5/C c h.) has a very small eigenvalue in 
the weak coupling limit. Hence, the solution of the integral equation (f4.52p 



\ V -p) = i r X _ Ar W , (4-55) 

' A-'cons. O/V/ch. 

is totally dominated by the small eigenvalue component. To see this, the unit operator 



13 To see that the second expression does produce the right speed of sound, one must know the explicit 
forms of I-p and I £ up to 0(XT 2 ). Since these forms are not essential to the present discussion, evaluation 



of the inhomogeneous terms I-p , I £ and the speed of sound v s are deferred to section IV E 
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l = £l/i></<l. (4-56) 

i 

where |/i)'s are the eigenvectors of (1— IC cons ,— 5/C c h.) with eigenvalues a*, may be inserted 
in Eq. ( |4.55| ) to yield 



l^> = E^I/*X/*W 

Cm 

= -|/ 5 ></ 5 |Xp>x(l + e>(A)). (4.57) 

In the last line, the eigenvector |/ 5 ) = |6 5 ) + 0(A) corresponds to the eigenvalue which 
vanishes in the A — > limit. The leading order value of a 5 can be obtained by 

«5 = (h\0-—^com-SK, ch .)\f 5 ) 

= -(b 5 \6JC ch .\b 5 )x{l + 0(X)). (4.58) 

As is shortly shown, due to the statistical factors, the eigenvalue a$ = 0(A) even though 
it originates from 0(X A ) correction to 0{\ 2 ) rungs. 

Using the leading order solution (|4.57| ), the leading order bulk viscosity is 

C = (3(lp\F\Vp) 

= -P (Ip ]f™t 5){ !? lIp) x (1 + 0(A)) • (4.59) 



At high temperature, the bulk viscosity explicitly contains mp hys provided by two factors 
of Xp. 

To obtain the explicit form of the ladder kernel including 0(X 4 ) rungs, note that when 
the external operators are bilinear, any diagram contributing to the Wightman function 
calculation can be regarded as a ladder diagram. The side rail part, as before, consists of 
two-particle intermediate states, and the rung part consists of the two-particle irreducible 
sub-diagrams between two side rails. Hence, if the rung matrix M. contains all possible 
rungs, and the propagators in T are the full propagators, the integral equation 

V p {k,q-k)=Xp{k,q-k) + J -^M(k-p)F(p,q-p)V-p(p,q-p) (4.60) 
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l 2 I3 I4 I5 I2 h h h 

FIG. 21. Diagrams corresponding to the C(A 2 ) 2-4 scattering processes. Mirror images and different 
permutations of l^s are not explicitly shown. Combining any two of these diagrams, including the square of 
a diagram, and integrating over U (i = 2, 3, 4) produce 0(A ) rungs. Due to the 4-momentum conservation, 



may correspond to the sum of all diagrams contributing to the correlation function of the 
bilinear part of the operator P. The exact solution of this equation, of course, is impossible 
to obtain. However, as shown in Eq. ( 4.59| ), the integral equation need not be solved; only 
the leading order number changing part of the kernel and the zero modes of the number 
conserving part are needed to evaluate the leading order bulk viscosity. 

To extract the relevant terms in the kernel, it is convenient to include in the "pinch- 
ing pole part" the rungs up to and including 0(X A ) corrections all calculated with free 
propagators 



/C pp = (Mo + M!)F pp . 



(4.61) 



Here, A4o includes only the one-loop rung, and A4\ includes C(A 4 ) rungs that can be 
obtained by adding two more lines to the diagrams in Fig. O, or equivalently, rungs that 



can be obtained by squaring the 2-4 amplitude shown in Fig. 21. M.\ also contains various 



corrections to the simple one-loop rung shown in Fig. |22|. Among these rungs, only those 
containing a number changing process are important in calculating the leading order bulk 
viscosity. The form of the pinching pole side rail matrix jF pp is the same as in Eq. ( [4.6| ) 
except that the self-energy Hj now includes contributions of 0(X 3 ) and 0(X A ) diagrams 



shown in Fig. With these definitions, the deviation from the "pinching pole part" arises 
only from the side rails. 
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FIG. 22. C(A 3 ) and C(A 4 ) corrections to the simple one-loop rung. The first two diagrams correspond 
to a part of the quartic vertex correction, the third one corresponds to the finite width correction. 




FIG. 23. C(A 3 ) and C(A 4 ) diagrams contributing to to the imaginary part of the A0 4 theory self-energy. 
Rewriting the integral equation as 

\Zp) = (l-lC ps ,-81Q\Vp) , (4.62) 
multiplying it with 5)C and adding the result to the original equation produce 

(l + 6JC)\Xp) = (I - !C PP - 5!C!C pp )\Vp) + 0((5!C) 2 ) . (4.63) 

Here 

8K /C pp = M 5F PP M T pp x (1 + 0(X 2 )) , (4.64) 

where <5jF pp is the C(A 2 ) correction to the free single particle spectral density in T pp . The 
non-pinching pole part does not contribute when sandwiched between two .Mo's as shown 
in appendix [FJ. 

The integral equation Eq. ( |4.63j ) is different from Eq. ( [4.52j ) in two ways. First, the 
inhomogeneous term in Eq. ( |4.63| ) is not purely Xp, and the kernel is not yet separated into 
the number conserving and the number changing part. Since only the leading order cal- 
culation is considered, the extra term in the inhomogeneous term is unimportant. For the 
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separation of the number conserving and the number changing part, it is more convenient 
to evaluate directly the needed matrix element 

(h\6IC ch .\b 5 ) = -(b 5 \(l - /C pp - 5/C/C pp )|6 5 ) , (4.65) 

rather than separating the two parts. 

Applying u T and using the same procedures as before, the integral equation (|4.63| ) is 
reduced to 

I' p {k) = D f {k) ~ j^h #buik& h) SUh) , (4.66) 

where 

I!p{k) = u T {k){l + 8K)Xp{k) , (4.67) 

and 

x ^ / 1[ SU-h)) ^(tu + k) (\%({k},k)\ 2 + T 3 ({k},k)) 

1 r 3 / d 4 l- \ 3 1 

+ j / II ( (2^4 S U-U) J Wdi^U + k) |T 4 ({/,},fc)| 2 J . (4.68) 

In the A0 4 case, the scattering amplitude involving 4 particles, Zi({7i},/c), includes the 
lowest order amplitude A, and 0(X 2 ) and 0(X 3 ) corrections. Since this part of the kernel 
conserves the particle number, the explicit form of 7l({li},k) is not important in the bulk 
viscosity calculation. 

The lowest order scattering amplitude involving 6 particles is given by 

%({h},k) = \ 2 Y,G% ee (k+lj+lk) , (4.69) 
where the sum is over all 10 different combinations of three momenta from the set {/,}. 
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(a) (b) 

FIG. 24. Diagrams corresponding to the squared terms in |7g| 2 . Figure (a) corresponds to the finite 
width correction the one-loop rung, and figure (b) corresponds to the finite width correction to the pinching 
pole part of the side rail. 

~ 2 

The square \T§\ 2 contains singular terms of the form G { ^ ee (li+lj+lk) which produces an ill- 
defined product of delta functions. However, these products of delta functions are removed 
by an additional term T 3 given by 

A 



T 3 {{k},k) = -Y.(Gl e \k+l 3 +l k )-GT(U+l 3 +h 



2 
A 4 

~Y 51 Pireeik+lj+lk) 



(4.70) 



The 0(X A ) rung diagrams corresponding to the terms in \%\ 2 can be obtained by attaching 
two more lines to the 0(X 2 ) rungs in Fig. [19] in all possible ways consistent with the A0 4 



theory. The prefactor 1/24 accounts for the symmetry factors of the diagrams. 

To see the connection between the terms in \%\ 2 and the diagrams, consider, for ex- 
ample, the 10 squared terms in \%({k} , fc)\ 2 + T 3 . Using the symmetry in l^s under the 
integral, these terms can be re-expressed as 



(4.71) 



T sq ({h},k) = 4A 4 [\Gf e (l 2 +l 3 +h)\ 2 + pfr ee (k+h+kY/2 

+ 6A 4 (\G% ee (h+l 2 +l 3 )\ 2 + Pil Uh+h+h) 2 /2 



The first term does not contain external momenta l\ and k. Hence, the cut lines correspond- 
ing to h,hyh form a cut self-energy diagram (together with the 4- momentum conserving 
5-function) corresponding to the finite width correction to the one loop rung shown in 
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Fig. |24|a. The symmetry factor associated with this diagram is 1/6. Factors in Eq. ( [4.68 ) 



combine to yield 1/12 including an extra factor of 1/2 from the relation of the form ([4.40 ) 



used to obtain Eq. ( [4.66] ). The remaining term containing h corresponds to the once 
iterated one-loop rung shown in Fig. |24|b. The symmetry factor associated with this di- 



agram is 1/4. Factors in Eq. Q4.68p combine to yield 1/8 again with the extra factor of 
1/2. 



Using the definition ( |4.22|) the above integral equation can be expressed symbolically 



as 



\l!p) = (1 - K hulk )\Dp) . (4.72) 

The self-energy £/ and the kernel -fTtmik has the following relationship 

/ ^J" n ( l °) S ^o(l)K^(l,p) = P^i(p) , (4.73) 
which can be proven by a similar argument used to obtain Eq. ( |4.49| ). From this relation, 



it is simple to see that the zero modes b^P) = / M E/(f) of the leading order ladder kernel 
(1 — K pp ) are still the zero modes of the modified ladder kernel (1 — K^ u \ k ). However, 
a previous zero mode corresponding to number conservation, 6 5 (f) = sgn(Z°) £/(/) is no 
longer a zero mode since 

J -j?L sgn(/°) n(l°) S free (/)iWZ,p) = sgn(p°) E 7 (p) + 5C(p) , (4.74) 

or symbolically, 

(6 5 |(l-#buik) = -(5C| , (4.75) 

where SC(p) contains only the number- changing part of the kernel. Only the overlap 
(5C\b 5 ) is needed to calculate the bulk viscosity since the leading order bulk viscosity is, 
ignoring corrections suppressed by O(X), 

( = P(Ip\F\V p )=P(I p \D p ) 

{fy\bs){bs\r p ) = jip\b 5 )(b 5 \i P ) 

P (h\(l-K hu , k )\b 5 ) ~ P (8C\h) ' 1 • ] 
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Straight forward calculations yield 



r d 3 \ 

(b 5 \Ip) = J [l+n(^)l <E l ) I P (E h 1) , (4.77) 



and 



-(5C\b 5 ) = (b 5 \(l-K hulk )\b 5 ) 
f 2 d% 

= 2 / II /n h ^12^3456 ^12 

x n(Et) n(E 2 ) [l+n(E 3 )\ [l+n(E 4 )] \l+n(E 5 )} [l+n(E 6 )\ , (4.78) 
where the differential scattering cross-section of 2 to 4 scatterings is given by |17j 

6 j3i / 

^12^3456 = g {27r) z 2E . (\ T ^ih})\ 2 ) ^YKli+U-U-U-k-h) I \mE 2 v 12 A\) , (4.79) 
where v\ 2 is the relative speed between li and — \ 2 



Vl2 



Jl + ii 

Ei E 2 



(4.80) 



In Eq. ( |4.79j ), all underlined momenta have positive energy. In the 2-4 scattering ampli- 
tude, products of delta functions do not appear since combinations of on-shell momenta 
in the propagators cannot be on-shell due to kinematic constraints. Hence, the T 3 term is 
irrelevant here. 

Using Fermi's Golden rule [17|, the expression — (5C\b 5 ) in Eq. ( f4.78|) may be interpreted 



as (2 times) the total 2-4 reaction rate per volume in a thermal medium. Due to the 
statistical factors, 0{m t ^ momenta dominate the integral in the expression ( |4.78|) . Hence, 
when T = 0(m phys /V\), (5C\b 5 ) = 0(A 3 T 4 ). Then since (b B \Ip) = 0(m 2 phys T 2 ) = 0(AT 4 ), 
the viscosity ( = 0(f3m phys / \ 3 ) = 0(T 3 /\) which is 0{\) smaller than the shear viscosity. 
Note that even at very high temperature, the bulk viscosity is a non-trivial function of 
the 0(1) ratio m 2 h /XT 2 . Hence, even at high temperature, the distinction between the 
thermal mass and the physical mass is important. 
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FIG. 25. Diagrams corresponding to the 0(A 3 / 2 ) 2-3 scattering processes. Mirror images and different 
permutations of ij's are not explicitly shown. Combining any two of these diagrams, including the square of 
a diagram, and integrating over k (i — 2, 3) produce C(A 3 ) rungs. Due to the 4-momentum conservation, 
I4 = h+h~h—h- 



D. Ladder summation for the bulk viscosity calculation with an additional gc^ 

interaction 

In the g0 3 +A0 4 case, the lowest order number-changing process is O(Xg) = 0(g 3 ) 
corresponding to 2 particles colliding to produce 3 particles. A few 0(A 3 ) rungs containing 
the effect of these scatterings are shown in Fig. [L3| Other rungs may be constructed 
by combining two diagrams among those in Fig. ^5] corresponding to the 2-3 scattering 
processes, or equivalently adding one more line to the rungs in Fig. IIU] in all possible ways 
consistent with the theory. 

With these rungs, the procedure used in the A0 4 case can be again used to produce the 
reduced integral equation (|4.66| ) now with 

^ bu i k aZi) = -^(l-e-^)(l-e-^) 

/ ft (tSi 5UHi)) (27r) 4 <KE /, + £;) \U{k}M 



2 



2 J i\\(2n 
(2tt) 



i=l 



(4.81] 
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The scattering amplitude involving 4 particles, ^({/j}, k), includes the lowest order ampli- 
tude A and O(Xg) and 0(X 2 ) corrections. Since this part of the kernel conserves number, 
the explicit form of %,{{h},k) is not important in evaluating the bulk viscosity. The lowest 
order scattering amplitude involving 5 particles is given by 

T 5 ({h},k) = Xg^G^ih+lj) - tg 3 J2G% ee (h+lj)G% ee (lm+Q , (4.82) 

where the first sum is over all 10 different combinations of two members from the set {k, k}, 
and the second sum is over 15 different combinations of four members of the same set. Since 
all the momenta in the set {h,k} are on-shell due to the delta function in Sf ree (— h), the 
combination k+lj cannot be on-shell. Hence, the ill-defined delta function products do 
not appear in \T 5 ({li}, k)\ 2 . Also for the same reason, it makes no difference whether the 
retarded propagators are used or the time-ordered propagators are used. For the sake of 
consistency, the retarded propagators are chosen here. The prefactor 1/6 accounts for the 
symmetry factors of the diagrams. 

The self-energy £/ and the kernel -Kbuik again has the following relationship 

/ (^ Pn ( /0 )^ e c(/)^buik(/,p) =£"£i(g) , (4.83) 

which can be proven by a similar argument used to obtain Eq. (|4.49| ). Again, the zero 
modes = L 2/(1) of the leading order ladder kernel (1 — K pp ) are still the zero modes 
of the modified ladder kernel (1 — -Kbuik) the remaining zero mode 6 5 (Z) = sgn(/°) £j(Z) 
is no longer a zero mode due to the number changing term in the kernel. The previous 



formula leading order bulk viscosity (|4.76 ) 



(SC\b 5) 

still holds with the same (b$\Ip) in Eq. ( |4.77 ) but with 



C = -P iIp l^ 5 \ Ip) , (4-84) 



(5C\b 5 ) = (b 5 \(l-K hulk )\b 5 ) 
1 r 2 d 3 l 

= 2J U^ d ^^ v ^n(E 1 )n(E 2 )[l+n(E 3 )}[l+n(E 4 )}[l+n(E 5 ^ . (4.85) 
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Here the differential scattering cross-section of 2 to 3 scatterings is given by 

dria-MB = II (27r) 32^ l r 5(fe})| 2 (27r)^a i+ ( 2 -(3-/ 4 -/ 5 ) / {^E^X) , (4.86) 

where t>i2 is the relative speed between li and — 1 2 . All underlined momenta in Eq. ([4.86 ) 
have positive energy. In the 2-3 scattering amplitude, the products of delta functions do 
not appear since combinations of momenta li+lj in the propagators cannot be on-shell due 
to kinematic constraints. 

Again, the expression — (5C\b 5 ) can be interpreted as (1/2 times) the total 2-3 re- 
action rate per volume. Due to the statistical factors, 0(m t h) momenta dominate the 
integral. Hence, when T = 0(m phys /\/X), (5C\b 5 ) = £>(A 5/2 T 4 ). The bulk viscosity is 
then ( = 0(T 3 /y/X) which is 0(X 3 ^ 2 ) smaller than the shear viscosity. Note that since 
the contribution of 0(m t h) momenta dominates the integral, the bulk viscosity is a non- 
trivial function of the dimensionless 0(1) ratio g 2 /Xm 2 h . Hence, even at high temperature, 
the distinction between the physical mass m p h ys and the thermal mass m t h is important. 
As shown in |1V E| below, the only place where m p h ys appears is in Ip. The propagators 
elsewhere must contain the thermal mass m t h- 



E. Inhomogeneous terms 

The explicit forms of the inhomogeneous terms and Ip(k) to leading order in 

weak coupling will be required in the following sections. For the shear viscosity, i^-(fc) 
corresponds to an insertion of the traceless stress tensor 7ii m given in Eq. ( [2.9|) . Hence, 

h{k) = hk m -\5 lm k 2 . (4.87) 

For the bulk viscosity, Ip(k) corresponds to an insertion of the operator V = V—v 2 e, and 
includes, as shown in section |111 D[ a one-loop "renormalization" contribution from chain 
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diagrams. Since the standard form of the stress-energy tensorj^] 



separates into the "kinetic" part, <9 M 0<9 V 0, and the Lagrangian part, g^C, the inhomoge- 
neous term for a pressure insertion may be separated into three parts, 

I P = 4 in - + I c + /£ hain , (4.89a) 

I e = J £ kin - -I c + J £ chain , (4.896) 

with the obvious notations. The sign difference in the Lagrangian term is due to the metric 
g^ = diag(-l, 1, 1, 1). 

The contribution from the "kinetic" part of V and e in the zero external momentum 
limit is simply 

= |k 2 , (4.90a) 

I^-(k) = E 2 k . (4.906) 

To determine the contribution of a Lagrangian insertion, Ic, it is convenient to rewrite the 
Lagrangian ( |1.1| ) using the equation of motion, as0 

£ = -EU] + -^-0 3 + , (4.91) 
2 m 2x3! 4L ' v ' 



where 



E[<j>] = <P^-£ = ( d2 - m l-^-^ 2 ) * . ( 4 - 9 2) 



Renormalization requires the counter term of the form ST^y — A^d^d^ — fj^u d 2 )4> 2 where A is a (infinite) 



constant |18|. However, this term does not concern us because its contribution to the inhomogcncous term 
is 5I^ V = (g^q 2 — <lp,<ly)A where is the external 4-momentum which is set to zero in the viscosity 
calculations. 

15 Alternatively, one can, of course, perform a straightforward diagrammatic analysis of the various terms 
arising from a Lagrangian insertion. 
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FIG. 26. The lowest order diagrams contributing to G r / > 2 C . Crossed circles indicate 4> 2 insertions, 
is times the equation of motion. An insertion of the operator E[<p] in a time-ordered 



iV-point correlation function simply produces an overall multiplicative factor [18 



N 

G(xi, ■ ■ ■ , x N ; E[4>(x)]) = K x ~ x a) G(x h ■ ■ ■ , x N ) , (4.93) 

a=l 

provided an irrelevant disconnected contribution is suitably subtracted (most simply by 
using dimensional regularization). To evaluate Ic, consider the correlation function of the 
Lagrangian L with some other bilinear operator, such as <fi 2 . In Euclidean space, 

G 02£ (x,y) = (T(0 2 (x)/:(y))) 

= ^(t(0 2 (x) E\m\)) + ^ T ^ x ) <^))> + 1 W 2 (^) ^y))) 



5(x- y ) {<p\x)) + ^ w>) nv))) + <P\y))) ■ 



(4.94) 



The lowest order diagrams for G^c are shown in Fig. [26|. The one-loop diagram in the figure 
is independent of external momentum. Hence, when the external frequency is analytically 
continued and the discontinuity taken, the contribution of this simple one-loop diagram is 
zero. The other two two-loop diagrams represent insertions of half the 0(XT 2 ) one-loop 
self-energy. Hence, to lowest order in the weak coupling limit, a Lagrangian insertion 
merely produces a vertex factor of 

I C = l^ml . (4.95) 
To determine the chain diagram part of the inhomogeneous term, evaluation of Re C-p(0) 



and ReC e (0) is required (c.f. Eq. (|3.14j)). The details of this evaluation are given in ap- 
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pendix |D|. The result, neglecting sub- leading corrections suppressed by 0{\f\), is simple, 



F**» = (A - 4") ReCp(O) = ~Sw^, (4.96a) 
m th 1 



/fain = (A _ jr\ ReCM = _ (4 _ 966) 



77? 



th 



Combining all three parts (and ignoring sub- leading corrections), the inhomogeneous 
terms for the pressure and energy density insertions are 

I v {k) = §k 2 , (4.97a) 

I £ (k) = E 2 k - 5m 2 th = k 2 + m 2 phys , (4.976) 

where 777 2 hys is the physical zero temperature mass. Note that to lowest order 777 2 hys is 
equivalent to m 2 h — (Tdm 2 h /dT)/2. 

Using these explicit forms of the inhomogeneous terms, the integrals in Eq. (|4.51|) can 
now be performed to show explicitly that the parameter v 2 in V = V—v 2 e is equal to the 
speed of sound v 2 = dV/de, 



v 2 = J (03^) [1+^(^)1 l k2 / / (^s^) \ l +<E k )\ (k 2 +mj hys ) 



1 5 m 2 



phys 



+ 0(\ 3/2 ) = v 2 s + C(A 3/2 ) . (4.98) 



3 12vr 2 T 2 

The details of evaluating the speed of sound, and the various integrals involved are given in 
appendix [D]. Note that even though the energy is defined with the thermal mass m tn , 
the speed of sound approaches 1/3 as the zero temperature mass m p h ys goes to zero, not 
as m tn goes to zero. In the massless limit, m p h ys — > 0, the stress-energy tensor is traceless 
due to scale invariance-p 6 ] In terms of equilibrium thermodynamic quantities this implies 
that e = 3V, and the speed of sound is f 2 = 1/3. 



16 Scale invariance is, of course, broken quantum mechanically, and this leads to a trace anomaly. This 
implies that the relation between the pressure and the energy density is modified to 3V— e — /3(A)T 4 /24 2 
where (3(X) — (TdX/dT) = 3A 2 /167r 2 [|l3) is the renormalization group /3-function. (The explicit forms of 
V and e are given in appendix [5]. ) Consequently, the speed of sound also receives an C(A 2 ) correction of 
8v 2 s = 5/?(A)/576tt 2 . 
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Combining the results of I-p, I £ and Vg, the inhomogeneous term for V insertion is 
(ignoring sub- leading corrections), 

Ip(k) = I v (k)-v 2 s I £ (k) 

/I 2\i 2 2 2 



where 



m 



phys f A (m th , T)k 2 - A 2 (m thj T) 



3 \A2{m th ,T) + ml hys A (m th ,T) 



, , 5k 2 



Khys {j^2T2 ~ 3 ) KW^ « 1) , (4.99) 



A n (m th ,T) = J -0^n(E p ) [l+n(E p )) |p|" . (4.100) 

Note that the inhomogeneous term is directly proportional to the zero temperature mass 
squared mp hys . Hence, Ip(k) vanishes in the massless limit, m p h ys = 0, and consequently, 
so does the leading order bulk viscosity (|4.21[) 



V. HYDRODYNAMICS AND THE BOLTZMANN EQUATION 

Kinetic theory and the Boltzmann equation have traditionally been used to calculate 
transport properties of dilute many-body systems when (except during brief collisions) the 
underlying particles can be treated as classical particles with well defined position, energy, 
and momentum. 



17 This is another consequence of (classical) scale invariance. Using the constitutive relation ( |2.l| ), the 
thermal average of the trace of stress-energy tensor may be expressed as (ignoring higher derivative terms) 

{Tp =W-e- CV-u , 

in the comoving (u(.t) = 0) frame. Classically, vanishes in the massless limit, and so docs its equilibrium 
average 3V—e. Hence, the leading order bulk viscosity £ must also vanish in the massless limit 
Again there are higher order corrections because scale invariance is broken quantum mechanically, and 
3V— e = 0(f3(X)). Consequently, in the massless limit, the bulk viscosity is ( = 0(XT 3 ). Note that 
when T 3> m p hy S /A, the massless limit estimate 0(\T 3 ) can be larger than the massive theory estimate 
0{<U>?T). 
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For this classical picture of particle propagation to be valid, the mean free path must 
be large compared to the Compton wavelength of the underlying particle. At extremely 
high temperature, the mean free path scales as 1/T. Hence, at high enough temperature, 
a Boltzmann equation describing the fundamental particles cannot be readily justified. 
However, as noted in previous sections, at high temperature the typical 0(l/\ 2 T) mean 
free path of thermal excitations is always larger than the 0{1/ \f\T) Compton wavelength 
of slowly varying thermal excitations. Consequently, as will be discussed, a Boltzmann 
equation description of effective thermal excitations with a temperature dependent thermal 
mass and thermal scattering cross section can be consistent at any temperature. 

However, there is a fundamental complication when attempting to formulate a Boltz- 
mann equation for effective excitations. In a non-equilibrium situation, the temperature 
may vary in space and time. Since the thermal mass depends on temperature, this implies 
that the mass of the effective excitations also varies in space and time. The correct treat- 
ment of this will be described shortly. As a warm up, first consider the usual Boltzmann 
equation with constant mass particles. It can be formulated from the statement that the 
rate of change in the comoving density of particles with an on-shell momentum k at posi- 
tion x equals the difference between the rates at which particles in this phase space region 
are generated or lost due to collisions, oiff] 

Here f{x,k) is the single particle density function, (df(x,k)/dt) gSLin is the rate of increase 
of the density of particles with momentum k at x due to collisions, and (df(x, k)/dt) loss 
is the corresponding rate at which particle density is lost. Here, _F ex represents whatever 
single particle external force or proper time derivative of 3-momentum, F ex = dk l /dr is 
present, and 7 is the usual E^/m. In this section an underlined momentum signifies an 



18 Here the momentum k is the canonical momentum, not the kinetic momentum. This is choice is 
necessary for the measure <i 3 x<i 3 k to be invariant along particle trajectories. 
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on-shell momentum with only positive energy. The collision terms on the right hand side 
may be expressed in terms of scattering cross sections and distribution functions. For 
instance, the gain of particles with k at x due to scattering is given by 

( df(x,k)\ 



I = E RfeJKlWU) fn(x;{p in })F m (x;{p o J,k). (5.2) 

gain in, out 



V at 

Here, T is a multi-particle scattering amplitude describing a process in which a particle 
with momentum k is produced in the final state. The initial probability density to find 
n-particles with momenta {p in } at position x is given by the product of single particle 
densities 

fn(x;{pJ) = Uf(^Pj (5-3) 

in 

under the assumption of molecular chaos. The factor F m (x; {p t },k) is the Bose or Fermi 
statistical factor for an m-particle final state with momenta {p out } and k, and is also given 
by a product of single particle statistical factors [l±f(x,p)] (with the upper sign for bosons, 
as considered here). Implicit in Eq. ( |5.2| ) is the assumption that the duration of scattering 
and the Compton wavelength of particles are short compared to the mean free time (and 
the scale of variation in the external force), so that the densities of particles participating 
in the scattering are accurately represented by the densities at a single position x. The 
loss rate {df(x,k)/dt) loss has a similar expression but with the particle with momentum k 
among the incoming particles. 

For simplicity, consider the A0 4 theory. If the interaction strength is weak, only the 
first few terms in the sum in Eq. (|5.2|) are important. The included terms must contain the 
leading order number- changing scattering processes since as discussed in earlier sections, 
the bulk viscosity calculation requires such scatterings. In that case, Eq. ( |5.1| ) combined 
with Eq. ( |5.2| ) becomes the basic relativistic Boltzmann equation, 



dt E k 7 '- / (2tt) 3 

r d 3 p 3 d 3 k 

~ / To \3 To \3 a<J l^3kVi2 
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(f(x, pj f(x, p 2 ) F(x, pj F(x, k) - F(x, pj F(x, p 2 ) f(x, p 3 ) f(x, jefj 
r d 3 pi d 3 p 2 

+ 2 / (n \3 In "^12-345*: Vn 

J12345 (2n) 6 (2n) 6 
(f(x, pj f(x, p 2 ) F(x, p 3 ) F(x, pj F(x, gg) F(ar, fc) 

pj p 2 ) /(:c, p 3 ) /(:c, pj f(x, pj /(:c, fc) ) 
/■ rf 3 p 5 rf 3 k 

+ / To VWo \3 acr 5fc^l234^5fc 

J12345 (ZTTj^ (2lT) d 

(ffapj) f(x,p 2 ) f(x,p 3 ) f(x,pj F(x,p 5 ) F(x,k) 

-F(x,p 1 )F(x,p a )F(x,^)F(x,p A )f(x,p s )f(x,$ > ) , (5.4) 



where, for a convenient presentation, the momentum space volume element <i 3 k/(27r) 3 is 
included in the formula, and a short hand notation F(x,p) = [l+f(x,p)} is used. The 
subscripts of the integral signs indicate that k is not integrated. Here, dcr 12 ^3fc is the usual 
two-body differential cross section, 

d 3 Pi 



da 



12^3fe 



II {2 n) 3 2E t ( 27r ) 4 %+g2-g 3 -^ ^(g^&.A) / (4E 3 E k v 12 2) , (5.5) 



where V\ 2 is the relative speed between particles with momenta p and — p , and the 
symmetry factor of 2 in the denominator arises from the non-distinguishability of the final 
particles. The 2-4 differential cross section da-12^3456 is given by Eq. ( |4.79| ). The factor 
of 2 difference in the second and the third term is due to the fact that the second term is 
symmetric in p 3 , p 4 , p 5 while the third term is symmetric in Pi, P2, P3? P4- 

The scattering interaction at a given point x still conserves energy and momentum 
even in the presence of a (slowly varying) external force which changes the 4-momentum 
of an excitation during the free flight between collisions. Hence, when multiplied by k? 
and integrated over k, the right hand side of the Boltzmann equation ( |5.4| ) vanishes. This 
implies for the left hand side 
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In the absence of any external force, this would become the local conservation equations 
for energy and momentum with the usual kinetic theory stress-energy tensor, 

With an external force, one instead finds 

d,TZ.{x) = S: x (x) , (5.8) 
where the source is given by 

wm-^™^- < 5 - 9) 

Now consider the case of excitations with a space-time dependent mass, m(x). Following 
the usual derivation of the equation of motion for a relativistic point particle, one finds 

F^(x) = j- k* = -d l m(x) , (5.10) 



where r is the proper time. Since the energy in this case, = yk 2 +m 2 (x), is space- 
time dependent, the partial derivatives in Eq. ( |5.6j ) cannot be simply taken outside the 
integral. Including additional space-time derivatives of energy changes the external source 
term ( p.9[ ) to 



= / (S 5 f{x ' m (it) ' m(x) / « r 



2^M)Jmmn**>- (5 - n) 



Finally, consider the situation of real interest where the mass of the effective excitations 
depends on the presence of other excitations in the medium. In the particular case of scalar 
interactions, relevant for A0 4 theory, the effective mass m(x) (to the lowest order) will be 

m 2 (x) = mj hys + 5m 2 (x) , (5.12) 

where 

^M/^W' (5 ' 13) 
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This is simply a rewriting of the one- loop result for the thermal mass m 2 h when f(x,k) 
is identified as the usual Bose factor. The g</> 3 +A0 4 theory result can be obtained by 
replacing A — ► X—g 2 /m 2 (x). From now on, 5m 2 (x) is assumed to have the form shown in 
( f5.13| ). In this case, one can define a modified energy- momentum tensor which satisfies the 
conservation equations ignoring 0(\ 2 ) corrections, 

T^(x) = J^—f( x ,® {p&+\gr8m\xj) . (5.14) 

At first sight, it may seem surprising that a conserved stress-energy tensor can be defined 
when the effective mass is space-time dependent. However, the underlying scalar theory 
does conserve energy and momentum. Hence, when the stress-energy tensor is correctly 
defined, the result must still be a conserved tensor regardless of whether the system is in 
or out of equilibrium. Appendix [D] shows that the form of T^ u in Eq. (|5.14|) is identical to 
the equilibrium expression of the field theory stress-energy tensor up to 0(A). Hence, the 



Boltzmann equation (|5T4]), with F* x (x) = —d l m t h(x), may be regarded as a kinetic theory 
description of effective temperature dependent excitations. 



A. The hydrodynamic limit of the Boltzmann equation 

Hydrodynamic excitations are arbitrarily long-lived, long-wavelength fluctuations that 
characterize near-equilibrium behavior of interacting fluids. Consider describing the re- 
laxation of a many-body system after a small initial disturbance. In a few mean free 
times, virtually all particles will have suffered numerous collisions with other particles in 
the medium. Hence, fluctuations in most degrees of freedom will relax in a few mean free 
times. However, an excess of a locally conserved quantity cannot simply disappear locally; 
to smooth a long- wavelength fluctuation in a conserved quantity, the conserved quantity 
must be physically transported over a distance comparable to the wavelength. For an 
arbitrarily large wavelength, this will require an arbitrarily long time. Consequently, for 
times long compared to the mean free time, the relaxation of the system may be described 
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solely in terms of very long-wavelength, long-lived fluctuations in locally conserved quan- 
tities. These are the hydrodynamic fluctuations of the system. For a simple fluid without 
an additional conserved charge, the only locally conserved quantities are the energy and 
momentum. 

To solve the Boltzmann equation ( |5.4| ) in this near-equilibrium hydrodynamic regime, 
the single particle density f{x,k) may be expressed as a small perturbation away from a 
"local equilibrium" distribution 

f(x,k) = f(°\x,k) + f {1 \x 1 k), (5.15) 

where the local equilibrium distribution f^°\x,k) is characterized by a local inverse tem- 
perature (3{x) and unit local 4-velocity u^(x) (satisfying = —1), 

f {0) (x,k) = 1/ (e-/^K(*)fc" _ i) . (5.16) 



To make the decomposition ( p. 15 ) unique, one must specify four conditions which serve 



to define the choice of local temperature and velocity. The most convenient choice is the 
Landau-Lifshitz condition |],|T!|, which fixes the local temperature f3{x) and the local 4- 
velocity u^(x) by requiring that f^°\x,k) produces the complete energy flow of the fluid 
|i"9| . This means that 



T^{x)u v {x) = T^{x)u v {x) 

= -e(x)u^(x) , (5.17) 

where T^ u (x) is the conserved kinetic theory stress-energy tensor which is defined in 



Eq. ( |5.14[ ), and 

d 3 k 



■ i ' [ i /. i i /■," i. ' — i r nil/ ' 

e(x)+V(x))u^(x)u u (x) + V{x)cT > (5-18) 



is the "local equilibrium" contribution to the stress-energy tensor, characterized by an 
energy density e(x) and pressure V(x). In other words, the correction f^'(x,k) to the 
distribution is required to make a vanishing contribution to T^ v u v . 
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By itself, the local equilibrium density function f^°'(x,k) is not a solution of the Boltz- 
mann equation. It makes the binary collision term in the Boltzmann equation ( |5.4j) vanish, 
while the convective derivative on the left hand side is non-zero (unless e and are con- 
stant). However, if e(x) and u^{x) vary on macroscopic scales, the size of the derivative 
|V/^(x, k)\/f^(x, k) will be small compared to any microscopic inverse length scale. 
Hence, by adding a correction to the local equilibrium distribution function, one may 
find a solution to the Boltzmann equation in which the size of the correction f^(x,k) is 
small, 

f {1) (x,k)\ ~ l ircc \Vf (0 \x,k)\ « |/ (0) (aU)| . (5.19) 

Here Zf ree ~ 1/na is the mean free path. 

To produce an equation for the first order correction f^'(x, k), it is convenient to choose 
the frame where u(x) = at some particular position x, so that 

f (0) (x,k) = n(x,E k ) = ePix) l k _ 1 , (5.20) 

and to express the first order correction as 

f (1 \x,k) = -n{x,E k ) [l+n{x,E k )]<f>{x,k) . (5.21) 

where (f)(x, k) is a slowly varying function to be determined by the (linearized) Boltzmann 
equation. Also, the lowest order temperature dependent part of the mass can be evaluated 

as 



5m 2 (x) = A- / - — — — n(x,Ek) 
1 



X—T 2 (x) (5.22) 
24 v ) \ j 

ignoring sub-leading terms. 

By equating the derivatives of the local equilibrium density function f^°\x,k) from the 
left hand side of Eq. ( |5.4| ) with the terms linear in <p(x,k) from the collision term, the 
following equation for <j)(x, k) in the u(x) = frame is obtained, 
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1 f dp 2 

= ^ y n 2^(2^)3 T ^PvP^l^^ (27i) 4 S(p l +p 2 ~p 3 ~k) 

x [l+n(z, ^)] [l+n(z,S 2 )]n(ar,S 3 )/[l+7i(a:,S fc )] 



+4/n d3p - 



MiP^k) (2vr) ^(Pi+Pa-Pg-^-Ps-i) 



24 7 ai 2^(2tt)3 
x Ei)] [l+n(x, J5 2 )] n(x, E 3 ) n(x, E 4 ) n(x, E 5 )/[l+n(x, E k )\ 

x (<f)(x,k) + 0(a;,p 5 ) + 0(x,p 4 ) + 0(x,p 3 ) - (j>(x,p 2 ) - <l>{x,pj) 



4 . 

T 6 ({p.},k) (2tt) +£2+23+24 -Pg-fc) 



48 7 ^ 2^(2^)3 

x [l+n{x,E x )} [l+n{x,E 2 )} [l+n(x,E 3 )\ [l+n(x,E 4 )]n(x,E 5 )/[l+n(x,E k )\ 
x {(f){x,k) + (j)(x,p 5 ) - 0(a;,2 4 ) - 0(^,Pg) - 0(^,P 2 ) - 0(^P X )) • 



(5.23) 



The inhomogeneous term on the left is a polynomial in momentum and derivatives of the 
flow velocity u(x), 



I(x,k) = P(x) (|k 2 -^ 2 (x)(k 2 +mj hys )) V-u(x) 

+ : ~7r^ f^%~|^i k2 ) (ViUj(x)+VjU i (x)-~5 ij V-u(x) 



(5.24) 

Note that choosing the u(x) = frame does not imply that the gradient at x, ViUj(x) 
are zero (but d^vP(x) = since = —1). Also note that I(x,k) contains mp hys not 



the thermal mass m 2 (x) 



m 



phys" 



-8m 2 (x). This agrees with the result of the previous 



section ( |4.99[ ) where the inhomogeneous term also lacked the thermal mass correction. In 
simplifying the left hand side (the inhomogeneous term), the equilibrium thermodynamic 
identity 



dT/T = dV/(e+V) 



(5.25) 



for the local thermodynamic quantities,P] and the lowest order energy-momentum conser- 



19 



This is a direct consequence of the form of the local distribution function /(°) ( 5.20| ) 
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vation equations 

d 

—e(x) = -(e(x) + V(x)) V-u(x) (5.26a) 

^n(x) = -VV{x)j(e{x) + V{x)) , (5.266) 

in combination with the definition v 2 = (dV/de) are used to rewrite the time derivatives 
of f3(x) and u(x) in terms of the spatial derivatives. 

Given the form of the inhomogeneous function I(x,k), rotational invariance (in the 
u(x) = frame) requires that <ft(x,p) have the following form, 

4>(x,p) = /5(x)vl(a;,p)V-u(x) 

+ (M— 1%) B ( X >P) (v l u J {x)+V J u l (x)-l8 ij V-u(x)) , (5-27) 

where p is the unit vector in the direction of p. Here, A(x,p) is the amplitude of the spin 
(divergence) perturbation, and B(x,p) is the spin 2 (shear) perturbation amplitude. The 
solution in any other frame is, of course, related to the given solution (j)(x, k) by a Lorentz 
boost. The linearized Boltzmann equation Eq. ( |5.23| ) is completely local in position; the 
parameter x is simply a label and henceforth will be omitted. 

The scalar process and the tensor processes decouple and can be studied separately. 
For the spin component, the integral equation for A is obtained by replacing I(x,k) in 
Eq. (pD by 

I P (k) = lk 2 -v 2 s (k 2 +m 2 phys ) , (5.28) 

and 4>(x,p) by A(p). For the the spin 2 component, Eq. ( |5.23| ) simplifies to the inhomoge- 
neous linear integral equation 

hk m -\5 ln y = i I n 2^(^r)3 ( 27 4 %+E 2 -23-^) \ r ^P v PA^)\ 2 

x [1+niE,)} [l+n(E 2 )]n(E 3 )/[l+n(E k )} 

x (B lm (k) + B lm (pz) - B lm (p2) - B lm {p^)) , (5.29) 

where Bi m (p) = (pip m —^Si m )B(p), since, as discussed earlier, the 2-4 scattering terms are 
unnecessary for the shear viscosity calculation. 
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After solving these linear equations to find the first order correction (p(x,p), the viscosi- 
ties can be evaluated by computing the first order correction to the stress-energy tensor 
and comparing it to the constitutive relation (|2.1|). One finds 



r r/ 3 k 



r r/ 3 k 

A(k)V-u{x) + - (kikj-iSij) B(k) (v i M j (x)+V j M i (x)-|5 ii V-u( 



and (using the lowest order result ViT^(x) = (s+V) ViU^(x)) obtains 



(5.30) 



77 = 15 I j2^fE k k ' n{Ek) ^ 1+n ^ B ^ ' ( 5 - 31fl ) 



r r/ 3 k 

C = PJ j^~ k (lk 2 +l5m 2 (x))n(E k )[l+n(E k )]A(k) 
r <7 3 k 

= N j2nfE k (l k2 -^(k 2 +m5 hys ))n(^,)[l+n(^)]A(fc) . 



(5.316) 



In the last expression, the Landau-Lifshitz condition Q5.17| ) Tm(s) = has been used to 
express the bulk viscosity in terms of the same source term that defines the amplitude 
A(k) (c.f. Eq. (pj). 



B. Equivalence of the Boltzmann equation to the field theory result 

The Boltzmann equation result for the shear viscosity calculation, Eq. ( |5.29|) , and the 
formulae from the finite temperature field theory Eq. ( |4.39| ) are remarkably similar. In 
fact, once the free cut propagators are used to force all momenta on-shell, Eq. (|5.29 ) will 



have exactly the same form as Eq. ( |4.39| ), provided that the thermal dispersion relation is 
used, and the scattering amplitude for temperature dependent effective excitations ( [4.47] ) 
is used in the collision term. 
With the identification of 
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D n (k) <-> {hk m -\5i m )^° ltz ik)B{k) (5.32) 

and 

'phys 



^ m phys + &m 2 , (5-33) 



the equation ( |5.29| ) for the spin 2 amplitude B(k) can be rewritten as 

J w = D n (k) - (1-e-^) / ^ WCfe-p) S'to.Cp) 2 ^\ p) > (5-34) 

where 5f ree (p) is the "free" phase space amplitude for an effective single particle thermal 
excitation, 

5 frcc (A;) = sgn(A; ) [l+n{k )} 1itb(k 2 + m t \) , (5.35) 

and the kernel Lsoitz. and the "self-energy" £p° ltz - are, 

L Bo it z .(k,p) = ^ J ^y^^ |^(p,P2;P3,&)| 2 (27r) 4 5(p+p 2 -P3-&) 

x S ircc (p 2 ) Sf tee (-p3) , (5.36) 

C'© = ^ (1-e-^) / -0- 4 L Boltz Xk-p)SUp) • (5-37) 

These are identical to the kernel L full (fc, p) and the self-energy Ej(/c) in the previous section. 
To write the equation in terms of a single function D w (x,p), an exchange of labels 3 <-> 1, 
and 2 ^> 1 has been used, together with the fact that the scattering amplitude squared 
| T4 1 2 is symmetric under the time-reversal. The only subtlety in deriving Eq. ( |5.34| ) is 
that it is written in terms of on-shell momenta with both positive and negative energies 
while the Boltzmann equation is written solely in terms of on-shell momenta with positive 
energies. The equivalence is possible because the negative energy contributions vanish due 
to the kinematic conditions enforced by energy conservation. 

With the use of the thermal mass in place of the zero temperature mass, Eq. ( |5.34 ) 
with the scattering amplitude 



T(PiP 2 ,P 3 :k) = A - g 2 (G R (p 2 +p) + G R (p 2 -p 3 ) + G R (p 2 -k)) , (5.38) 
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has exactly the same structure as Eq. (|4.39|) . 



For the bulk viscosity, the equivalence can be shown by first noting that the right hand 
side of the integral equation ( p.23|) also has a near zero mode corresponding to letting 
(f)(x,p) be a constant, or integrating over d 3 kn(Ek)[l+n(Ek)}/ (2n) 3 Ek. In particular, if 
the identification 

Dp(k) <-> Zf oltz ik)A{k) , (5.39) 



is made, 0(x,p)'s in the right hand side of Eq. ( |5.23|) are replaced by 2 (to ac- 



count for both the positive and the negative energy on-shell momenta), and integrated 
over d s kn(Ej ; .)[l+n(Ei i: )]/(27T) 3 Ek, the result on the right hand side is identical to 
(6 5 |(l-^buik)|& 5 ) = ~(SC\b 5 ) in Eq. ( j!7g ). The integral equation ( g2g ) itself is not 



strictly identical to Eq. (fl66|). Eq. ( |4.66|) includes additional pieces in the number con- 
serving parts of the kernel. However, the number- changing part of the two kernels in the 
equations ( f4.66| ) and ( p.23| ) are the same. 



The existence of the near-zero mode implies that the integral equation for the spin 
component A is again dominated by this near- zero mode component. By using the same 
arguments as in section [IV Q the leading order bulk viscosity can be again written in terms 



of the near-zero mode matrix element. Since Ip(k) is an even function of the definition 
of the inner product of two functions ( |4.22|) can be used to express the bulk viscosity as 



r d 3 k 

C = PJ j^- k Ip(t) n{E k ) [l+n(E k )] A(k) 



where (b 5 \Ip) and (5C\b 5 ) are again given by Eq. (|4.77|) and Eq. ([4.78|) . 

As stated earlier, the Boltzmann equation for the fundamental particles ceases to be 
valid at temperatures high enough to make the mean free path smaller than the Compton 
wavelength of the underlying particle. Thus, it may be surprising to find that a form 
of Boltzmann equation remains valid at all temperatures, provided the mass parameter 
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is interpreted as the mass of effective thermal excitation and a temperature dependent 
scattering amplitude is used. 

It is interesting to consider the effect of using the thermal mass and the thermal scat- 
tering amplitude in various temperature range. At low temperatures (T<^m p h ys ), the 
thermal correction to the mass (c.f. Eq.( |5.13] )) is completely negligible. Consequently, for 



the leading order calculation, thermal quantities are unnecessary and the viscosities may 
be calculated by using kinetic theory of the non-relativistic particles. If the temperature 
is in the range m p h ys < T^m p h ys / V% most particles are highly energetic, but the ther- 
mal corrections to the mass and the scattering amplitude are negligible. Consequently, 
the viscosities at these temperatures can be calculated by the kinetic theory of relativistic 
particles with the zero temperature mass and the scattering amplitude. At very high tem- 
peratures, T3>m p h ys /A, all mass scales, including the cubic coupling constant, other than 
temperature may be ignored. Hence, both the shear and bulk viscosity may be calculated 
from the kinetic theory of massless excitations with only the quartic interaction.^ 

The most interesting region is at intermediate temperatures T = 0(m phys /y/\). Since 
this is much larger than m p h ys , and the typical size of loop momenta at high temperature 
is 0(T), one might expect that the replacement of zero temperature mass by the thermal 
mass should have a negligible effect. This is true for some observables, such as the shear 
viscosity. However, for the bulk viscosity the contribution from momenta of (9(m t h) is not 
negligible compared to the hard momentum contribution. 

To understand the behavior of the bulk viscosity, first note that the classical scale 
invariance requires the classical bulk viscosity to be proportional to m phys /T ||. When 
T = (9(m p hy S /v / A), the effect of quantum mechanically broken scale invariance is negligible 
compared to the m p h ys . In a scalar g0 3 +A0 4 theory, the 2-3 amplitude for soft momenta 
is T 5 ~ 0(Xg/m^ h ) Hence, the expression (SC\b 5 ) above is a non-trivial function of the 



20 The power counting performed in this paper is also valid for the massless scalar theory since the 
excitation in this case develops non-zero thermal mass of 0(y/\T) at non-zero temperature. 
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dimensionless 0(1) ratio g 2 /Xm 2 h . As discussed in section |TV D| , (6C\b 5 ) = 0(A 5 / 2 T 4 ) in 
this temperature range. Consequently, 

C = ^ d h ^{g 2 /\m 2 th ) x (1 + 0(m th /T)) , (5.41) 

where c?b u ik is a dimensionless function of 0(1). The coefficient function rfbuik cannot be 
calculated from the massless scalar theory. Thus, including the thermal correction to the 
mass is essential to calculate the correct leading weak coupling behavior of bulk viscosity 
when T = 0(m p hy S /\/A)- 

If the temperature is high enough, T ^> m p h ys / A, the quantum scale anomaly dominates 
the effect of the physical mass term, and the leading weak coupling behavior of the bulk 
viscosity is identical to that of the massless theory with only the quartic interaction. Due 
to the scale anomaly, (b^\Ip)(h) is non-zero but proportional to the 0(X 2 ) /3-function, 
(h\I P )(k) = C(/3(A)k 2 ). Then since (5C\b 5 ) = C(A 3 T 4 ), C = 0(AT 3 ) when T > m phys /A. 
Note that at this temperature, C(AT 3 ) is larger than (9(m 4 hys /A 3 T). 

In contrast, the shear viscosity is insensitive to the soft momentum contribution. For 
the typical momentum k = 0(T), g 2 G(k) = 0(g 2 /T 2 ) < 0(X 2 ). Hence, in this case, the 
scattering amplitude is dominated by the quartic interaction term, or, T ~ A. Then, the 
dimensional analysis demands that 



T 3 

V = rfshear , (5.42) 



where c4h ea r is a pure number of 0(1) which can be calculated from the kinetic theory of 
massless excitations with only the quartic interaction. 



VI. CALCULATION OF VISCOSITIES 

The purpose of this section is to apply the results of previous sections to the calculation 
of the bulk and the shear viscosities in the A0 4 and the g<j) 3 +\(j) 4 theory. We begin with 
the bulk viscosity. Recapping the result (|4.T6|) , the leading order bulk viscosity is given by 
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(5C\b. 

Among the factors in Eq. ( |6.1| ), 

d 3 l 



C = -P ih ti ib A Ip) ■ (6-1) 



f d i 

(h\Ip) = J —^[l+n{E l )]n{E l ) I P {E h \) (6.2) 

reduces to a 1-d integral over the magnitude of the momentum, and can be easily evaluated 
numerically. The denominator for the A0 4 theory 



- .M. 

12^3456 V12 



i (2vr) s 

x n{E x ) n(E 2 ) [l+n(E 3 )] [l+n(E 4 )} [l+n(E 5 )] [l+n(E 6 )} , (6.3) 
with the differential cross section 

^12^3456 = II TTxkw l^6({I i })| 2 (27r) 4 5(Z 1 +Z 2 -Z 3 -Z 4 -Z 5 -Z 6 ) / (AE.E^AV) (6.4) 

i= 3 i 2 ^r 2E i I 

has a complicated angle dependence through the scattering cross section Tq. Using the 
energy-momentum conserving delta function and rotational invariance, 4 of the 18 dimen- 
sional integral involved in calculating (5C\b 5 ) can be done. The remaining 14 integrals 
involve 5 integrations over the magnitudes of momentum, and 9 angle integrations. Due to 
rotational invariance, one of the solid angle integration can be trivially done reducing the 
expression to a 12-dim integral which must be evaluated numerically. The denominator 
for the g0 3 +A0 4 theory is similarly given but with the more complicated 2-3 scattering 
amplitude (c.f. Eq. (|4.85| )). In this case, a 9-dim integral must be carried out numerically. 



In contrast to the bulk viscosity calculation, the shear viscosity calculation requires 
solving an integral equation. Here, the integral equations derived in section fV| for the 
effective vertices are further reduced to one dimensional integral equations, and the explicit 
form of the A0 4 theory kernel is evaluated and briefly examined. The shear viscosity in 
terms of the effective vertex D n is given by 

V = TR i 7^1 n(k°) S iTCC (k) , (6.5) 



10 J (2rr) 4 wv ' v y ee v ' Ej(fc) 
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where D n satisfies 

I n (k) = D^k) - (1-e-^) / ^ LUt,p) SUP) • (6-6) 

The kernel Lf u n is 

1 f dH d A l 2 
L fn a{k,p) = -J Jc^i SirUh) S ilee (-l 2 ) (27r) 4 5 (h-h+p-k) T{ly,p,l 2 ,k) , (6.7) 

and involves the g0 3 +A0 4 tree level scattering amplitude 

T(k,P}k,k) = A - g 2 (Gnh+p) + G R (k-k) + G R (k-l_ 2 )) . (6.8) 

The imaginary part of the self-energy is 

Mt) = ~(l-e-^°) / L iun (k-p)SUp) , (6-9) 

where the free cut propagator contains the thermal mass, 

Sfree(p) = [l+n(p )] sgn(p°) (2rc) 4 -5(p 2 ! +m t 2 h ) . (6.10) 

The inhomogeneous term I n (k) = k[k m — |5; m k 2 represents an insertion of the traceless 
stress tensor 

lTi m = di<f>dm<p ~ \5lmd k <P k (j) . (6.11) 

Eq. ( |6.6| ) is a set of 5 independent (due to the traceless symmetric spatial indices) 
3-dimensional (since all momenta are on-shell) linear integral equations. Using spatial 
rotational invariance, this may be further reduced to a single one dimensional linear integral 
equation. In section |IV|, the effective vertex D^ijp) was shown to be an even function of 
the momentum. Isotropy then requires that D n (p) have the structure 

Ar(p) = (piP m -h8im)D sheaI (\p\) , (6.12) 



where p is the unit vector in direction of p. Contracting Eq. ( |6.6| ) with kik m and carrying 
out the frequency integration with the help of the on-shell 5-function yields a single one- 
dimensional integral equation for -D shear (|p|) 
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n /ii i\ i 2 , 3 f d \p\ P 2 

£>shear(|k|) = k 1 



■n(E p ) [l+n(E p ))N sheai (\k\, |p|) ZW(|p|) , (6.13) 



16 J (27r)»JSgr p 

where the kernel N shear is an angular average of the 3d kernel Lfuii, 

iV shear (|k|,|p|) = (1-e"^) J #dcosfl(cos 2 fl-§) 
x ~L Ml (E k , k; E p , p) (e^-1) - L full (£ fe , k; -E p , p) (l-e~^) 



(6.14) 



Here, # is the angle between the vectors k and p. As before, the on-shell thermal width 



T p is defined as 



£/(£ P ,p) 
2E n 



(6.15) 



Note that iV shear is symmetric, A^hear(|k|, |p|) = iV shear (|p|, |k|), since L Ml (E k , k; E p , p) = 
e /3(B fe - Bp ) L{ull (E p ,p;E k ,k). 

Finally, in terms of the scalar function Z} s hear(|k|), the shear viscosity is given by the 
simple integral 







60vr 2 



d\k\ 



n(E k )[l+n(E k ))D sheai (\k\) 



(6.16) 



The final one-dimensional integral equation (|6.13|) must be solved numerically. Clearly, 
one must first evaluate the full g0 3 +A0 4 theory kernels N s h eai ( |6.14| ). With both cubic and 
quartic interactions, the full "rung" L fu n (p.7|), is too complicated to compute analytically. 
However, in a pure A</> 4 theory, the integral is reasonably straightforward and one finds, 
4 A 2 c dH 

Lf ull (k-p) = — J -j—-^ S bee {-1) Sireeil+k-p) 



A 2 [l+n(A; u -p )] 
87r/3|k-p| 

x I 6((k-p) 2 ) In 



1 — exp(— (3r + {k— p)) 



+ 0(-(k-pY-Am 2 th ) In 



1 — exp(— (3r_(k— p)) 

sinh ( f3r + (k—p) / 2 ) 



sinh(/5r_(fc— p)/2) 



(6.17) 



where 



r ± (k-p) = - (|k-p|^/l+4m t 2 h /(^-p) 2 ± (k°-p°] 



(6.18) 
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The term involving 9 ((k—p) 2 ) in ( |6.17|) represents the result of integrating over two on- 



shell delta functions with Bose factors n(Ei)[l+n(Ei + i c -p)]. These statistical factors indi- 
cate that this term describes the transfer of an incoming momentum k—p to a thermal 
excitation with momentum [, producing another on-shell excitation with 4-momentum 
(Ei + k- p , 1+k— p) with a stimulated emission factor of [l+n(Ei + k- p )]. 

The second term involving 9(— (k— p) 2 — 4m^ h ) describes the usual process of creating 
two propagating on-shell particles with total invariant mass larger than twice the mass of 
the initial single particle excitation. It differs from the zero temperature result, 



\im LUk-p) = X 2 e(k°~P°) 9(-(k-p) 2 -Am 2 phys )^l+Aml hys /(k-p) 2 /lQ7r , (6.19) 

only because of the stimulated emission in the final state. 

To solve the integral equation ( |6.13| ) numerically, the magnitudes of momenta |k| and 
|p| need to be discretized in order to turn the integral into a finite number of coupled 

,4 

linear equations. Given the explicit form of the A0 4 theory kernel Lf ull (k— p), evaluating 
the coefficients of the linear equations requires numerically computing two 1-dimensional 
angle integrations; one for the angular averaged kernel A?^ loar (|k| , |p|), and the other for 
the self-energy Ej (p) (c.f. Eq. (|6.9|)). In contrast, for the g0 3 +A0 4 theory calculation, 
evaluating each coefficient of the final matrix equation requires first computing two 2- 
dimensional angular integrations since Lf u n(k,p) is no longer just a function of k—p due to 
the non-trivial structure of the scattering amplitude. 

Although the angular integrations involved in Nf hear is too complicated to carry out 
analytically, some qualitative behaviors of the kernels can be easily found. One property 
of the kernels, important in carrying out numerical analysis, is that Af s h ear (|k|, |p|) has 
a discontinuous first derivative across |k| = |p|. To see this, consider, for example, the 



expression of A^ ear (|k|, |p|) in Eq. ( p,14| ). If the angle 9 is defined to be the angle between 



two vectors k and p, the azimuthal angle integration in Eq. ( |6.14| ) is trivial. For the 



cos9 integration, it is convenient to change variable to y = |k— p|. The Jacobian of this 
variable change cancels an explicit l/|k— p| factor contained in . Then Eq. (|6.14j) may 
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be rewritten as 

Al 4 ear (|k|,|p|)= r^\dyh{E kl E p) y)-F + {E k ,E p ,y)) , (6.20) 
J ll k l-|p|l \ / 

where F_(E k , E p ,y) contains the logarithmic part multiplying 9((k—p) 2 ) in Lf uYl with 

(cos 2 9— ~), and F + (E k , E p , y) contains the logarithmic part multiplying #(— 4m 2 h — (fc— p) 2 ) 

in Lf ull with (cos 2 #— ~). 

Differentiating with respect to |k| yields, 

— jV s f ear (jk|, |p|) = (F_(E k , E p , |k| + |p|) - F + (£ fe , £ p , |k| + |p|) J 

-sgn(|k|Hp|)(F_(£7 fcj E p J|k|-|p||)-F + (E fc ,E p ,||k|-|p||)) 
+ J^ dy^(F4E k ,E p ,y)-F4E k ,E p ,y)^ . (6.21) 

Note that in the second line the signature of |k| — |p| explicitly appears as a result of 
differentiating the lower-limit of the integral. Since F±(E k , E p , | |k| — |p| |) is non-zero in the 
|k| — > |p| limit, this implies that the first derivative of A^ ear ( |lc| , |p|) has a discontinuity 
across |k| = |p|. 



VII. NUMERICAL RESULTS 



A. Choices of parameters and discretization method 



In this section, the 1-d integral equation ( |6.13| ) reproduced here for convenience, 



A hcar (|k|) = k 2 + ^/ M ^ v n{E p ) [l+n(E p )\ N Bheai (\k\, |p|) D shsai (\p\) , (7.1) 

is numerically solved to obtain the shear viscosity in the A0 4 theory, and the integrals 
involved in the bulk viscosity calculation 

[pC\h) 

are numerically carried out. 
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dx 
d|p| 




10 15 20 25 30 35 40 45 50 

IpI 

FIG. 27. A typical form of the Jacobian |dx/d|p||. Here, the parameters are M = 50T, fi = T. 



First consider the shear viscosity calculation. By discretizing the magnitude of momenta, 
the integral equation can be turned into a finite set of linear equations which can be 
straightforwardly solved by computer. 

The discretization method chosen is the lowest order two-point Newton-Cotes formula 

H, 



f XN+1 dx f(x) = |(/(x0 + f(x N )) + £ fl?i) + 0((Ax) 

Jx A i=2 



(7.3) 



where Ax is the distance between two data points. The reason behind choosing this simple 
discretization is following: Due to the discontinuity in first derivatives of the kernel (a 
kink), second derivatives at |p| = |k| is not well defined. Hence, in choosing a discretiza- 
tion method, higher order formulae are not necessarily more useful than the lowest order 
formula. 

To successfully implement the numerical analysis, a suitable parameterization must be 
chosen so that the improper integrals in the integral equations become proper ones. The 
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parameterization used here is basically a logarithm of a Fermi distribution with chemical 
potential M 

_ l n ( e (M-/(|p|))/3 + ]j 

X ~ \n(e M ? + 1) ( } 

where 

/(IpI) = IpI-t4x- + ca*- (7-5) 
IpI 

The new parameter x varies from 1 to as |p| varies from to oo. Constants c, M, and /i 
are adjustable parameters. 



The parameterization ([774 ) is chosen for the following two reasons: (a) To account for 
the soft-momentum contributions there must be enough data points near |p| = 0. (b) 
The kink in the kernel implies that the contribution from the momenta |p|~|k| cannot be 
ignored even for large |k|. This implies that a discretization that sparsely samples large 
momentum values are not suitable since it will miss the non-negligible wiggles in those 
regions. One must choose a parameterization that distributes data points more or less 
evenly in momentum space until the cut-off is reached. As shown in Fig. |2?| the Jacobian 
(which may be interpreted as the density of sampled points when Ax is constant) 

dx 

^l + e (/(|p|)-M)/*)ln(l + e J 



l + C/i 2 /(| P |+/i) 2 



d\p\ 

shows that these two conditions are satisfied by the parameterization given in Eq. (|7.4j ). 
The positive parameter c controls the height at |p| = 0, and thus controls the percentage 
of data points sampled below |p| = fi, where \i is usually chosen to be 0(m t h)- In the 



present calculations, c is chosen as 

2(M - 10 n) 



(M > 10/i) . (7.7) 



9/i 

This particular choice of c puts about ten percent of the total number of data points below 

IpI = n- 

Once |p| exceeds fi, the Jacobian is almost constant until |p| = M is reached. This 
implies that the data points are evenly distributed throughout the momentum range 

85 
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FIG. 28. Numerical solutions for D shcar (|k|) for m th /T = 0.05,0.5, 1.5. 

fi < |p| < M. Since the Jacobian drops off very sharply once the momentum exceeds M, 
one can control the largest momentum sampled by choosing the value of M. Usually the 
largest momentum sampled is the size of M. Any M much larger than any of the mass 
scales in the integrand will do. For the present calculation, the values of M ranges from 
30T to 100T as the mass increases. 



m = 0.05T 

m = 0.5T 




B. Numerical results 

For the shear viscosity, the value of the thermal mass m t h are chosen to be 0.01T < 
m th — 3T to show the high temperature behavior of the shear viscosity. Typical forms 



of solutions for .D s h ea r(|k|) are shown in Fig. |28| for m t h/T = 0.05,0.5, 1.5. The values of 



viscosity, extracted from the Newton- Cotes formula 

Voo = Vn + a/N 2 , (7.8) 
with the number of data points N = 100, 200, 400, 600, are shown in Fig. EUl As expected, 



the shear viscosity rises as the mass is increased. In the non-relativistic limit of large 
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rj\ 2 /T 3 
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FIG. 29. Upper curve represents numerical results for shear viscosity. Lower curve represents one- loop 
calculation result. 
The fit for the upper curve is rj(m t h/T) = 2860(T 3 /A 2 )(1.0 + 0.596m th /T + 0.310m t 2 h /T 2 ). And the 
fit for the lower curve is rj (m th /T) = 733(T 3 /A 2 )(1.0 + 1.33m th /T + 0.627m t 2 h /T 2 ). 



m P hys/T, as discussed in section [II B| , the shear viscosity must rise as 0( 



[m 



5/2 , 
physJ 



In the 



ultra-relativistic or high temperature limit, the viscosity is nearly independent of the mass 
and 0(T 3 ). In the high temperature limit, the result of resummation is about 4 times 
larger than the one- loop result alone. 

For the bulk viscosity, the numerical evaluation of the multiple integral involved in 
calculating (5C\b§) was carried out by Monte-Carlo method. Typically, each point in the 
plot Fig. |3D| is evaluated by about one million data points. As expected, the bulk viscosity 
behaves like 1/A 3 for small values of m t h/T and rises sharply as the ratio increases. 

VIII. SUMMARY 



Hydrodynamic transport coefficients can be evaluated from first principles in a weakly 
coupled scalar field theory at arbitrary temperature. Using the diagrammatic rules derived 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



m th /T 

FIG. 30. The plot of the numerical results for the bulk viscosities from resummed ladder diagrams. 
The fit is \og((\ 4 T/m* hys ) = 12.9 + 1.43(m th /T) + 1.36 log(m th /T). The error bar is typically about 1% 
of the the value of £. 

in ||, it was shown that an infinite number of diagrams contribute to the leading weak 
coupling behavior of the viscosities. The dominant diagrams were identified by counting 
the powers of coupling constants, including those generated by near "on-shell" singularities 
cut-off by the single particle thermal lifetime. An infinite class of cut "ladder" diagrams 
were found to make the leading order contributions. The geometric series of cut ladder 
diagrams was summed by introducing a set of effective vertices, satisfying coupled linear 
integral equations. These equations were reduced to a single integral equation, which 
was then shown to be identical to the corresponding result obtained from a linearized 
Boltzmann equation describing effective thermal excitations with temperature dependent 
masses and scattering amplitudes. The effective Boltzmann equation is valid even at very 
high temperature where the thermal lifetime and mean free path are short compared to 
the Compton wavelength of the underlying fundamental particles. 

Spatial isotropy allows one to reduce the dimension of the resulting integral equations 
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FIG. 31. The plot of the numerical results for the bulk viscosities from resummed ladder diagrams for 
m th < T. The fit is (\ 4 T/m 4 phys = 1.4x 10 6 ((m th /T) 2 + 0.18(m th /T) 4 + 0.12(m th /r) 6 ). Since m 2 jT 2 ~ A, 
the bulk viscosity at high temperature behaves like mp hys /A 3 T as expected. The error bar is typically 
about 1% of the the value of (. 

to one-dimension, at which point they must be solved numerically. Numerical results for 
the viscosities in a scalar A0 4 theory is reported. 
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APPENDIX A: THERMAL PROPAGATORS 

The imaginary-time single particle propagator is 

G E {k,iu n ) = 1 . — -, (Al) 

ul + k 2 + mf h + S E (k, iu n ) 

where S B (k, iu n ) is the full Euclidean self-energy, uj n is the discrete frequency 2nnT, and 
mth is the thermal mass containing 0(XT 2 ) thermal corrections. The Euclidean propagator 
has the spectral representation, 
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^ /, • x f dou p(k, uj) 

J 2,71 UJ — lOJ n 

Hence, the single particle spectral density is obtained by analytically continuing G E in 
frequency and taking the discontinuity across the real axis, 

p(k) = -i (G E {k, k°+ie) - G E (k, k°-ie 

—i i 
+ 



k 2 +m t 2 h + E(A;) k 2 + m t 2 t + T,(ky 

2E/(fc) 
|P + m t 2 h + £(A;)| 2 ' 



(A3) 



where 



S E (k, A;°+*e) = = £«(£;) - iZj(k) , (A4) 

is the analytically continued Euclidean self-energy. 

The real-time propagator {T(4>(x) 0(0))) used in the cutting rules is 

C Ws /|[l + „ MM | 4a )( M; ^_ it) j . (A5) 

By changing the integration variable uj to — uj and adding the two expressions together, 
the real-time propagator can be re-expressed as 

= -1 (G B (k, A;°+*e) + G E (k, fc°-ie)) + i coth(/3A;72) 

- - ^ (A6) 



A; 2 + m t 2 h + £(A;) A; 2 + m 2 h + £(£;)* ' 

where the fact that the spectral density p{k) is an odd function of the frequency is repeat 
edly used. 

APPENDIX B: EXPLICIT FORM OF THE LADDER KERNELS 



To analyze the zero momentum, small frequency limit of the viscosities, a detailed 
understanding of the structure of the ladder kernel /C is needed. The 4x4 kernel is a 
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product of two factor, JC = M.T . First, consider the A0 4 theory kernel. The 4x4 rung 
matrix is 



M(k—p) 



1 -iC(k-p) 

iC(k-p)* 
L(p-k) 

L(k-p) 

which has entries consisting of the uncut rung 



V 



(Bl) 



dH 



and the cut rung 



\2 

C(k-p) = -%— . 
V 11 2 J (2tt) 4 



L(k- P ) = ^ J *L S(l+k-p) S(-l) 



G(l+k-p)G(l) 



:b2) 



(2tt) 



(B3) 



As before, (5(7) is the uncut propagator defined in Eq. ( |A6| ) and S(l) = [l+n(/°)]p(7) is the 
cut propagator. 

The side rail factor of the ladder kernel is given by 



T{p, q-p) 



■ (B4) 



/ G(P) G(q-p) S(-p) S(p-q) G(p) S(p-q) S(-p) G(q-p) ^ 

S(p)S(q-p) G(-p)*G(p-qy S(p)G(p-q)* G(-p)*S(q-p) 

G(p)S(q-p) S(-p)G(p-qy G(p)G(p-qy S(-p)S(q-p) 

S(p)G(q-p) G(- P yS(p-q) S(p)S(p-q) G{-p)* G{q-p) 

When the external momentum vanishes, the matrix — p) can be written as a sum 
of four outer products 



F(P, -P) = w{p)u T {p) + h{p) f{p) + «(p)f T (p) + n(p)( T (p) 



(B5) 



where 



w T (p) = (l, 1, (l+e-^)/2, (l+e p0/3 )/2)[l+n(p°)]n(p ) p(p) 



MP) 



u 



\p) = (1, 1, (l+e^)/2, (l+e-v°P)/2) 



(B6) 
(B7) 



91 



h T {p) = (0, 0, 1/4, -e^/4) 4iL , (B8) 

f{p) = (0, 0, 1, -e-^) , (B9) 

K T (p) ee (1, e"^, e"^, 1) [l+n(p°)] 2 /b 2 + ml hys + S(p)] 2 , (BIO) 

?(p) ee (-1, -e"^, -1, -e"^) , (Bll) 

ee (1, e^, 1, e^) n(p°Y/[p 2 + m 2 hys + E^)*] 2 , (B12) 

( T (p) ee (-1, -e^, -e^, -1) . (B13) 

In section [IV] it is asserted that the hj T part is orthogonal to the inhomogeneous 
terms, and if T = T—hj T , then j T M.T = TM.h = 0. Showing that and are 
orthogonal to the inhomogeneous terms Ia(p) and z^(p) is trivial because both h(p) and 
have vanishing first and second elements, while the only non-zero elements of 1a{p) 
and za (p) are the first and the second ones, respectively. To show that hj T is orthogonal 
to M.T and TM., first note that h(p) is orthogonal to u(p), and ((p), and that 
j(p) is orthogonal to w(p), k{j>) and Hence, JF/i = j 7 T = 0. Due to the relation 

L(k—p) = e^ k °~ p0 ^L(p—k), h(p) and j T (p) are "eigenvectors" of the rung matrix, 



M(k-p)h(p) = (0,0, L(p-k), -e p0l3 L{k-p)) T '■ 



p{p) 



4£j(p) 

l0.0.L-e^) r L( P -fc)- p(/M 



4Ej(p) 

oc /i(jfe) (B14) 



and 



f{k)M{k-p) = (0,0, L{p-k),-e~ k0/3 L{k-p)) 
= (0,0,1, -e-^)L(p-k) 

oc /(p) . (B15) 

Hence, fMf = TMh = 0. 

In section IV, the relation (|4.15|) for the A0 4 theory rung, repeated here 



K pp (k,p) = u T (k)M{k-p)w(p) 



l-e- k P) L(k-p) S iicc (p) j 2£/(p) , (B16) 
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was important in simplifying the expression for the pinching pole contribution. This rela- 
tion can be proven as follows. Let w(p) = (1, 1, (l+e~ p0/3 )/2, (l+e p0/3 )/2). Then, 

u 1 (k)M{k-p)w{p) = 2lmC(k-p) 

+ (1 + e ^)(l + e" fc0/3 ) (L{k-p) + e-^- k °^L(p-k)) U 
= (l-e- k °P)L(k-p)(e p °P-l)/2 , (B17) 

where to obtain the last expression, the optical theorem 

Im C(k-p) = -(L(k-p) + L(p-k))/2 , (B18) 



and a symmetry of the cut rung, 

A 2 r d A l 

L{p-k) = — J — — S bee (l+p-k) S bee (-l) 



e^-^ p L(k-p) , (B19) 



are used. The second expression in Eq. ( |B19j ) is a consequence of the property of a cut 



propagator Sf ree (— k) = e~ k ° ^ S{ ree {k) . When combined with the remaining factors forming 
w(p), this yields Eq. ( |B16|) . Note that the real part of C(k—p) makes no contribution to 
the pinching pole part of the rung matrix. 

For the g<p 3 rung matrix representing the straight single line rungs, exactly the same 
argument applies to yield 

K Vmc (k,p) = u T (k)M Yme (k-p)w(p) 

= g 2 (i-e-^) S(k-p) S lrec (p) 1 2E x (p) , (B20) 

where the rung matrix A4n ne (k— p) is now given by replacing — iC with —g 2 G, and L with 
g 2 S. And also in place of the optical theorem for the imaginary part of C, the straight 
rungs satisfy 

ReG(k) = (S(k) + S(—k))/2 . (B21) 
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V/ 



FIG. 32. Diagrammatic representation of the rung matrix A^f u ii-box(^— p)- 



For the analogous relation for the box sub-diagram rung ( [4.37[ ), reproduced in Eq. ( |B23[ ) 
below, it is convenient to consider the "rung" matrix generated by (full) box sub-diagram 

dH 



M 



full-box 



(k-p) 



(2tt)< 



■ Mm*(k-l)r(l, -l)M Vmc (l-p) 



(B22) 



illustrated in Fig. The box sub-diagram rung MiH(k-p) is the non-pinching pole 
contribution of ■M{u\\-hox(^~P) ■ The separation of the pinching pole contribution and non- 
pinching pole contribution can be made after contracting with u and w. Since the line 
sub-diagram rung matrix is diagonal, contracting with u{k) and w{p) produces 

dH 



u (k)M{ U n-b ox (k-p)w(p) 



^(k)M^ e (k-l)^{l,-l)M^>(l-p)w^ . (B23) 



To simplify this expression, first, note that the elements of the rung matrix satisfy, 



M^_ hox (k,p) = e-^MS hoK (k,p) , 
M^_ hox (k,p) = e^M^ hox (kp), 



M^_ hox (k,p) 
M^_ hox (k lP ) 



full-box 



ik,p) 



^° + ^M^_ hox (k,p) 



(B24a) 
(B24b) 
(B24c) 
(B24d) 
(B24e) 
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(B24f) 



due to the corresponding property of the cut propagator S(—p) = e~ p °^S{p). Using the 
facts that the imaginary part of a uncut diagram is given by the sum of all possible cut 
diagrams divided by a factor of (—2) (c.f. Eq. (|3.11|) ), one finds 



u T (k) .Mfuii.boxfe P)w(p) 



1-e 



(B25) 



Note that since the frequency integral is not yet carried out, this relation is valid for both 
the pinching pole contribution and for the non-pinching pole contribution., 

For the full g0 3 +A0 4 theory rung matrix, the relations ( B24 ) can be again shown to 
hold. Hence, so does the relation (|4.40 ) 



u T (k)M {ull (k,p)w(p) = (l-e-^) L M1 (k,p)S free ( P )/2Z I (p) 



(B26) 



where 



WM = {M^(k,p) - e k0 ?M^(k,p) 



(B27) 



APPENDIX C: ZERO MODES OF LADDER KERNELS 

The integral operator (1— /C) where /C = has four zero modes as uo — > corre- 

sponding to the momentum-energy conservation. To see this, one must know the contri- 
bution of an momentum-energy density T 0fl insertion. The standard stress-energy tensor 
is given by 

= d^d^ + g^C , (CI) 

where the Lagrangian is given in Eq. ( |1 . 1| ) . The momentum density T 0i = d°<f)d' l <f) contains 
only the "kinetic" part, and hence its contribution in the zero spatial momentum limit is 
simply 
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j(k, u-k) = (0, F(£; -^/2), 0, 0) , (C2) 

where k is the loop momentum flowing through the lines connected to T 0t and u is the 
external frequency. 

The energy density on the other hand contains both the kinetic part d°(j)d (j) and the La- 
grangian part. The contribution from the kinetic part of T 00 in the zero spatial momentum 
limit is again simple, 

z^(k,u-k) = (0,k°(k°-Lu/2),0,0) . (C3) 

The contribution of an Lagrangian insertion can be calculated by applying the method 
used in the main text (c.f. Eq. flj.91|) ). The result up to 0{\ 2 ) is 

zl(k,u-k) = -(0, E(k-uj)*, 0, i[l+n(w-k°)}J: I (w-k)) (C4) 

where 

E(jfe)* = ^([l+n{k )}i:{ky -n(k )Z(k)) , (C5) 

is the uncut two-loop self-energy. 

Given the form of z % and z° = z^ in +z^, a tedious but straightforward calculation yields, 

V^{k,uj-k) = z*(k,u-k)T{k,u-k) 

= g*{k,u-k)/u- J -^g^p,u-p)M(p-k)T(k,u;-k)/u 

= g»(l-K)/u, (C6) 

where, the row vector is given by 

g"(k, u-k) = fc"(0, i{G{k-u)* - G{k)*), -iS(k),iS(u-k)) . (C7) 

Since z^T is finite as to goes to zero, 

lim g"(l -K) = 0. (C8) 
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Hence, in the zero external momentum, and zero frequency limit, the operator (1 — JC) has 
four left zero modes given by g^.f^\ Note that the first element of is always zero for 
any uj. Hence, trivially, <? m Xa = for all to. This implies that the inhomogeneous term is 
orthogonal to the zero modes of the operator lim(l — /C). 

Since V M corresponds to an insertion of the energy-momentum tensor, these zero mode 
can be used to verify that correlation functions of the energy or momentum density vanish 
in the zero spatial momentum limit. Full correlation functions involving must vanish 
as the momentum goes to zero since the conservation equation relates the time derivative 
of a conserved "charge" to the divergence of its current. Hence, for example, the correla- 
tion function of two momentum densities must behave like k 2 in small momentum, finite 
frequency limit. In terms of the effective vertex, this implies that 

V^V A = , (C9) 

for an arbitrary (non-singular) external operator A. Since, the effective vertex D a can be 
expressed as 

Va = ihc lA ' (C10) 

for momentum density correlation functions, 

= q^Xa/u = , (Cll) 

since as explained earlier q^Za = 0. 

For the pinching pole part of the integral equation, 



21 The right zero modes of (1 — !FM) / M can be also obtained in an entirely similar way. They are 
f"{k,w-k) = k"(i{G{k)-G{w-k)),0,-iS(w-k),iS{k)) . 



97 



I B (k) = D B {k) - (1-e-^) J ^ L {ull (k-P) SUP) , (C12) 

the (left) zero modes are 

b ll (k) = k ll [l+n(k°)]S ilcc (-k) . (C13) 

To verify this, note that [l+n(/c )] = 1/(1— e~ fc0/3 ) cancels the prefactor (1— e~ fc0/3 ), and 

/d^k . . 1/" d^k d^l\ d^li , , j , , , im2 

j^yL K Siree(-k) L^k-p) = -J ^-^ (2tt) 8(h+p-l 2 -k) \T(h, p; l 2 , k) \ 

x Sf ree (—k) Sf Tee (li) Sf Tee (—l 2 ) 

= 2p li E I (p)n(p°), (C14) 

where to obtain the second line, the original expression is averaged with two equivalent 
expressions differing by the k <-> or k <-> / 2 labeling changes. The scattering amplitude 
is, as before, 

Tih^k) = (A - <? 2 (GRih+p) + G fl (Zi-A0 + Gfl^-Za))) . (C15) 

Hence, when b^(k) is applied to the right hand side of the above integral equation, the two 
terms in the right hand side cancel each other exactly. 

APPENDIX D: STRESS-ENERGY TENSORS AND THE SPEED OF SOUND 

In this appendix, the equilibrium A0 4 theory stress-energy tensor is calculated, including 
the C(AT 4 ) correction. When the cubic interaction term is added to the Lagrangian, to 
this order, one only has to make the change 

q 2 

A^A-Ar. (Dl) 

In equilibrium, the stress-energy tensor is diagonal in the comoving frame, 

(THe q = diag( £ ,7>,P,P). (D2) 
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Then due to the equilibrium thermodynamic identity 



dT \T 

only the pressure needs to be calculated 



e = T 2 ^-(l-V) , (D3) 



FIG. 33. The lowest order connected vacuum diagrams for the thermal correction to the pressure in 



the g4> 3 +X(f> 4 theory. The cross indicate an insertion of 8m 2 a2 



The easiest way to calculate the correction to the free-particle pressure is to sum the 
contribution of all connected vacuum graphs (see, for example, Ref. 0]). To express the 
pressure in terms of the thermal mass, the Lagrangian may be rewritten as 

-c = i#-a» - v 2 + m t \)0 + A 4 _ j (D4) 

where m 2 h = m^+5m 2 h . To the lowest order, the zero temperature mass mo may be 
identified as the physical mass m p h ys - Then the first order correction to the free-particle 
pressure arises from the diagrams in Fig. where the term ^5m 2 h (f) 2 in the Lagrangian 
is treated as an additional interaction term. In Euclidean space, the free-particle pressure 
through 0(A) is, 

d 3 k k 2 1 



(2tt)3 3 ool + V + ml 



(2tt) 3 ujI + k 2 + m 



th, 



■ trn 2 th r d 3 k 1 

2 (2 7 r)3< + k 2 + m 2 1 ' 1 " J 



where u n is the usual discrete Euclidean frequency. 



Using standard techniques (see, for example, Ref. fl22[), the sum over u n can be sepa- 



rated into the zero temperature contribution and the non-zero temperature contribution. 
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At high temperature, the (renormalized) 0(mp hys ) zero temperature contribution is negli- 
gible compared to the 0(XT A ) thermal correction. Hence, to calculate the pressure through 
0(XT A ), only the non-zero temperature contribution needs to be examined. For the one- 



loop diagrams in Fig. [J^, the non-zero temperature contribution is obtained by a simple 
replacement of 

^ 1 n(E k ) , s 

n=0 U n + h k 

where n(Ek) is the usual Bose factor, and the energy here is El = k 2 +m 2 h . The thermal 
part of the pressure is then, 

Pth = / <E k ) (±k 2 + \8ml) , (D7) 



where the fact that 



A r <7 3 k 

*< = 2 (D8) 

to lowest order is used to simplify the expression. 

To evaluate the pressure ( |JJ7| ) , consider first the the one- loop thermal correction to 
the mass ( [L>8| ) . Since 5m% h involves an explicit factor of A, only the leading order term is 
needed. By setting the mass in the integrand to zero, the leading order contribution to the 
integral in Eq. (P8|) can be evaluated as 



/ d\k\ |k| Y e 



MP 



4tt 2 



E 4 = ^ , (D9) 



4vr 2 £J s 2 24 



neglecting sub- leading terms suppressed by 0(V\). 
Next, consider the integral 

d 3 k , „ s 1 



p °= 1 ww k n{ - E ^¥' '■ (D10) 
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The integrand depends only on the magnitude of the loop momentum |k|. Changing the 
integration variable |k| to m t hsinh# yields 



07T Jo 



xdrnx.it cosh 9 



m 



4 oo 
th 



E 



2tt 2 ^ V s2m t\ 



T 2 1 

K 2 {s(5m th ) + -K (s^th) 



(Dll) 



where i^„(x) is the modified Bessel function of order n. expression in Eq. (|D11| ) is obtained 



by integrating by part, and using the standard expression for the Bessel functions [23). The 
leading and the next-to-leading terms of Vq in the small parameter /3m t h can be calculated 
by using the small x expansion of K 2 (x), 

2 1 



x- 



0(x lnx) 



(D12) 



Substituting Eq. ( pi2| ) into Eq. ( [01 1| ) and performing the elementary sums yield, 

vr 2 T 4 m 2 T 2 



90 



6 th J 

24 



+ 0(T<) • 



(D13) 



All together, the thermal pressure up to C(AT 4 ) is 



th 



vr 2 T 4 m t 2 h T 2 T 2 ^ AT 2 



90 
tt 2 T 4 
90 



24 

2 T 2 
^phys- 1 

~24 



48 V 24 J 
T 2 (XT 2 " 



48 V 24 



(D14) 



From this, the thermal energy density (|U3[) can be easily calculated, 



£th = 3P t h + 



m phys^ 2 



12 



■ n 



(E k ) (E% - \5m 



(D15) 



To obtain the last expression, Eqs. ( |D7| ) and ( |D8|) are used. From the expression of the 
pressures ( p7| ) and the energy density (P15|) , the thermal expectation of the stress-energy 
tensor up to 0(A) can be compactly written as 

d 3 k 



eq 



(2tt) 3 Eu 



n 



(D16) 
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From the expressions ( pi4| ) and ( pi5| ), the speed of sound v" 2 = dV/de can now be 
straightforwardly calculated up to 0(m 2 hys /T 2 ), 



, _ (dV/dT) _ 1 5m 2 hys 



s (de/dT) 3 12tt 2 T 2 1 ' 



APPENDIX E: NEAR SOFT AND COLLINEAR SINGULARITIES 



In this appendix, the soft momentum behavior of diagrams contributing to the cal- 



culation of the viscosities is briefly examined. In section |TJ it is asserted that the near 
soft singularities do not affect the power counting described in that section. Here a brief 
demonstration is presented. The temperature is assumed to satisfy T ^> m p h ys so that 
ml = 0(XT 2 ). 

To the leading order, the effect of adding one more rung M.T is to provide the integrand 
one more factor of 

u T (k)MUk-p)w(p) = (l-e- k °?) LUk-p) [l+n(p )] S ircc (p) j 2E x (p) , (El) 
together with an additional integration over the 4- momentum p. Here Lf u \\(k— p) consists of 



the cut diagrams such as those in Fig. [19]. In section |T|, this additional factor is regarded 
as 0(1) since the inverse powers of A from the inverse of the self-energy is canceled by 
the explicit 0(X 2 ) scattering amplitude squared contained in Lf u \\. When the (on-shell) 
momenta k and p are soft, this argument could be upset if (a) the size of the self-energy 
is smaller than C(A 2 T 2 ), (b) the would-be soft singularities (factors of 0(T/m t h)) 
from the Bose factors are not compensated by the small momentum space volume, or (c) 
the non-pinching pole contribution of the side rail matrix is comparable in size to the 
pinching pole contribution. 

To see that none of these possibilities actually occur, first consider the size of the thermal 
"scattering amplitude" which is contained in the expressions for S/(p) and Lf u n(k, p), 



T(k,P}l*Q = X ~ 9 2 {Gr(Li+P) + Gnik-k) + G R {l x -l_ 2 )) . (E2) 
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As before, the arguments of the propagators in Eq. (|E2|) are all combinations of two on-shell 
momenta. The 4-momentum squared of the sum of two on-shell momenta satisfies 



(k±p) 2 - (E k +E p 



= 2[E k E p t k-p + m 2 th ) 
>2(E k E p -\k\\p\-m 2 h )+4m 2 th 



for all k and p since 

E k E p -\k\\p\-m 2 h = ((£ fc £ p ) 2 -(|k||p|+m t 2 h ) 2 ) / {E k E p +\k\ \p\+m 2 th ) 
= m t 2 h (|k|-|p|) 2 /(E k E p +\k\\p\+m 2 h ) 
> . 

Similarly, the 4-momentum squared of the difference of two on-shell momenta is 



(E3) 



(E4) 



(k±p) 2 - (E k -E p f = 2{E k E p =f k-p — m 



th) 



> 2(E k E p - |k||p| -m 2 th ) 
>0, 



(E5) 



for all k and p. Hence each propagator in (|Ei|) is bounded by l/m 2 h for all on-shell 
momenta [ { , k and p, so that 



\T(k,&kM = 0(X) + 0(g 2 /m( h ) = 0{\) , 

since by assumption g 2 = 0(\m 2 hys ) and (^ph ys / m th) — 1- 

With this estimate of the size of |T| 2 , the size of the g0 3 +A0 4 theory cut rung, 



(E6) 



Lm(k,p) = 7T 



1 r d% d% 



Stree(h) S iree (-l 2 ) (27i) 4 5(h-l 2 + P -k) T(l u p;l 2 ,k) , (E7) 



2 J (2vr) 4 (2tt) 4 

at soft k and p can be determined. When the external momenta k and p are 0(mth), the 
two 5-functions in the cut propagators can be satisfied by 0(mth) loop momenta. Hence, 
the two momentum integrations over these two (^-functions and the energy-momentum 
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conserving ^-function are 0(1). Consequently, the size of L{ n u(k,p) at soft k and p is 
determined by the size of the scattering amplitude squared |T| 2 = C(A 2 ) and two 0{T / m th ) 
factors from the statistical factors. Hence, all combined, 

Wm th ) = 0(\ 2 T 2 /m 2 th ) = 0(A) , (E8) 

since m t \ = 0(AT 2 ). 

Given this result, one may also estimate the imaginary part of the self-energy 

= \(l-e-^) J -0^ S {ree (k) Lwife, k) , (E9) 

at soft external momenta p. When k = 0(mth), the momentum integration together with 
^-function contribute a factor of 0(m 2 h ). The 0(T/m t h) Bose factor and the prefactor 
(1— e~- >/3 ) = 0(m th /T) combined are 0(1). Since L full (p, k) at soft k and p is 0(A), 
putting all terms together yields, for soft p, 

E f (p) = 0(Am t 2 h ) = 0(A 2 T 2 ) . (E10) 

Hence, £/(p) is C(A 2 T 2 ) for both hard and soft p. 

With all the ingredients at hand, the size of the soft momentum contribution to the 
additional rung u T (k)M.{ u \\{k— p)w(p) can be readily examined. When all the momenta 
involved are soft, the momentum integration combined with the 5-function in the cut 
propagator provides a factor of 0(m 2 h ). Once again, the prefactor (1— e~^ 0/3 ) combined 
with the Bose factor in the cut propagator is 0(1). The cut rung is L lua {k,p) = C(A), and 
the self-energy remains 0(X 2 T 2 ). Hence, all combined, the integration over u T M.w can 
be regarded as C(Am 2 h /A 2 T 2 ) = 0(1). Consequently, one can conclude that the power 
counting performed in section |ITJ is not altered by soft momentum contributions. 

For the near collinear singularities, note that the estimate for the scattering amplitude 
( |E6| ) holds for all on-shell momenta. Hence, there is no large factor resulting from near 
collinear singularities and the power counting performed in section [TIT] is again not altered. 
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APPENDIX F: DETAILS OF CHAIN DIAGRAM SUMMATION 





FIG. 34. Diagrammatic representation of the correlation function <jaa- Filled bubbles represent the 
sum of all cut chain diagrams. The external operator A is represented by black circles at each end. 

For the sake of simplicity, A0 4 diagrams are examined first. Diagrammatically, the 
sum of all A0 4 chain diagrams for the Wightman function a AA of a bilinear operator A 



can be represented by Fig. |34[ The filled bubble in Fig. [3J with the ends of the cut lines 
on opposite sides represents sum of all cut chain diagrams with the same topology. This 
sum is denoted by L cha _ in (q). Similarly, the filled bubble with the ends of the cut on the 
same side represents the sum of all cut diagrams with the equivalent topology. This sum 
is denoted by C c hain(g) 
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FIG. 35. Diagrammatic representation of the equations satisfied by the resummed cut and uncut 
bubbles. In this diagram, the external operators at both ends are (f> 2 rather than A. The sum of all chain 
diagrams for correlation function o aa can be expressed in terms of the solutions of these equations. The 
empty cut bubble is denoted by L(q), the empty uncut bubble C(q). The filled bubble on the left hand 
side of the first diagram is denoted by L c hain(g), the filled bubble on the right hand side of the second 
diagram is C ch ain(g)- 



The equations for L chain (q) and C^^q) are 



L chain (q) 






( 

+ x 




, Cchain(<?) , 




K C (q) J 


V 



C (qY L Q (q) \ ( L chain (q) ^ 



J 



Cchain(O') 



(Fl) 



Lo(-q) C (q) 

as illustrated in Fig. The cut bubble L (q) and the uncut bubble C (q) are again 
defined by Eq. (g^) and Eq. (^) with I A (l, q-l) = 1. 

The above matrix equation is easy to solve. The solution is 

Lo(q) 



J chain 



(q) 



and 



C c hnm(q) 



(l-AReC(g)) 2 + (f) 2 (L(g)-L(-g))2 
C{q)-\{{ReC{q)f-\{L{q)-L{-q)f 



AReC(g)) 5 



mm 



L{~q)f 



(F2a) 



(F2b) 



where the finite temperature optical theorem 



1 



Im C{q) = -- (L(q) + L(-q)) 



(F3) 



is used to obtain the form shown in ([F3). 

When cubic interactions are included, the "chain" diagrams also include cut two-loop 
diagrams shown in Fig. |9] where the bubble in the diagram may be regarded as the sum of all 
A0 4 chain diagrams. Equivalently, the vertex contribution may be regarded as containing 
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I A (l, -/) = I Ail, -I) + ARe CU(O) x (l + 0(VX)) . (F4) 

A straight forward application of cutting rules then yields, 

*S 9 (0) = -4^(0) (5(0) C A (0) + 4^(0)* (5(0)* L A (0) 
+ 4(^(0)* 5(0) C A (0) + AL A (0) 5(0) L A (0) 

2 4 

= _ 8 \ReC A (0)L A (0)+A^ T (Re C A (0)) 2 L o (0) , (F5) 



where the optical theorem Q3.ll]) , and 5(0) = g 2 Lq(0) / mf^, justified below, are used. 

For the remainder of this appendix, the estimates ReCU(0) = 0(T 2 ), ReCo(0) = 
0(1/ V\), 9 2 G(0) = 0(g 2 /m 2 th ), and 5(0) = C(1/A 2 T 2 ) used in this appendix and sec- 
tion |IIID| are examined. To estimate Re C A (0), consider the following explicit form of the 
real part of an uncut bubble C A (0) 

ReC A (0) = ~J ^Ll A (l,-l)Re(^G 2 (l)) 



dH 



I A (l, -I) coth(/°/?/2) 



47 (27r) 4 ' V ^[/2 + m 2 h + S (/)]2 [/2 + m 2 h + S(Z) *^ 

id r d 4 l 



I A {1, -I) coth(|/°|/?/2) 2tc5{1 2 + m 2 



4 dm 2 J (2tt) 4 
where to obtain the last line, the elementary relation 

d ( 1 \ 1 



x(l + 0(A 2 )), (F6) 



dx \x + a J (x + a) 2 
and the single particle spectral density in the weak coupling limit, 

P(l) 



(FT) 



[Z 2 + m t 2 h + £(/)] [P + < + E(l)*] 
= Pfrcc (/) x (1 + 0(A 2 )) , (F8) 

are used. Note that the identification m 2 = m 2 h must be made after the derivative is taken. 
Again at high temperature, the zero temperature contribution is smaller than the thermal 
contribution. Hence, the coth(/3^/2) factor in the integrand can be replaced by the Bose 
factor 2n(Ei) to calculate the leading weak coupling behavior. 
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For the pressure and the energy density insertions, 



1 d 
ReCV(0) = - 



d 3 l 



2 dm 2 J (2tt) 3 Ei 
1 d 



2 dm 2 



n{E l ) ±1 2 



T 2 



m' 2 =m?, 
th 



48 



(F9) 



and 



ReC £ (0) 



1 d 



d 3 \ 



2 dm 2 J (2tt) 3 E 1 



n{Ek) El 



th 



1 a I m 2 T 2 ' 
3V <> + — 



l2 = m ?h 



48 ' 



(F10) 



2 dm 2 ' 

ignoring higher order contributions. 

The expression for ReCo(O) is also given by Eq. (|3.8|) by setting Ia = 1- At high tem- 
perature, the leading 0(T 2 ) contribution to the integral in the last expression in Eq. ( p.8|) 
comes from loop momenta of 0(T). However, when the derivative with respect to m 2 is 
taken, this is zero. The next largest contribution to the integral is 0{mT) coming from 
loop momenta of 0{m). When the mass derivative is taken and the identification m — > m t h 
made, this yields ReC (0) = 0{T/m th ) = 0(1/ VX). 

To estimate the size of propagators at zero momentum requires knowledge of the size of 
the self-energy at zero external momentum. The lowest order imaginary part of the off-shell 
self-energy comes from the gcj) 3 one-loop diagram shown in Fig. 0. To estimate the size 
of this diagram, the propagators in the loop must be regarded as resummed propagators. 
Then the self-energy in the zero 4-momentum limit satisfies 



lim lim 2{l+n(q )] E/(g) = g 2 L (0) x (1 + 0(A)) 

= 0( 9 y\ 2 ) = 0(ml hys /X) 



(Fll) 



where pinching pole approximation of L(0) is again used. The above estimate for the 
self-energy implies that, 

lim lim s(9) = lim lim 
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^#x(l + (VA)) 

m th 

(% 2 /A 4 T 4 ) = 0(mJ hys /A 3 T 4 ) 



0(1/A 3 T 2 x (g 2 /T 2 )) < £>(1/A 2 T 



2rj-i2\ 



(F12) 



and 



lim lim g 2 Re G(q) = lim lim 



g 2 [q 2 + m 2 th + E R (q)] 



*0 q^O 



?°-o q-o |g 2 + m? h + E(g)| 2 
0(^ 2 /m t 2 h ) = 0(AmJ hys /m t 2 h ) < 0(A) 



(F13) 



APPENDIX G: THE IMAGINARY PART OF TWO-LOOP SELF-ENERGY 



Using the cutting rules, the imaginary part of the two-loop self-energy can be expressed 



as 



Sr" l00P (£) = ^(l-e-^)%), (Gl) 



where 

r dH d 4 k 

«%) = J J^Xi J^J SbeeO) S {vee (k+q) S iree (-l~k) . (G2) 

Using the on-shell ^-functions free particle cut propagators 

Sfree(0= E (Tiil+nfaEftirSQ - aM) / E t , (G3) 

o-i=±l 7 



frequency integrations in Eq. ( |G2|) can be straightforwardly carried out to yield 

1 r d 3 ld 3 ~k 

= a k+ g (Ti+k(Ti W7 —r E / — — [l+r2(a fe+9 E fc+g )] [l+n{aiEi)} 

x [l+n(a k+ iE k+ i)} 8(E q -aiEi-a k+q E k+q -a k+ iE k+ i) . (G4) 



The argument of remaining ^-function can be satisfied only when two of the a are +1 and 
the other one is —1. By suitably changing labels, the above then becomes 
3 r d 3 l d 3 k 

S ^ = / tp tp f ^ 1+n ^ [l+^l-^fc+g)] n(E k+ i) 5(E q +E k+ i—E k+q —Ei) . 

o[Ztt)° J hih k+q h k+ i 

(G5) 
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To carry out remaining integrations, the angles between spatial vectors are defined as 

cos6 k = k-q/|k||q| (G6) 
cos0j EEk-l/|k||l| . (G7) 

Changing variables to E k+q and E k+ i with the Jacobians 

dE k+q /d cos 6 k = \k\\q\/E k+q (G8) 
dS fc+I /dcos0, = |k||l|/S fc+ ,, (G9) 

one finds 

-%) = o/o^ 5 | i [dE l( l\k\ [l+n(Ei)} j Ekl dE k+l f Kq dE k+q [l+n{E k+q )\ 
s{Zir) |q| J JE kl JE kq 

x n(E k+l ) SiEq+Ew-E^-Et) , (G10) 

where 



£±EE^/(|k|±|l|)2 + m t 2 h , (Gil) 

with analogous definitions for E^ q . Carrying E k+q and E k+ i integrations amounts to figuring 
out the kinematic conditions. Straightforward calculation then yields 

f pE q / e -PEi -/3(E q +E k )\ 

x | dEk [l+niE^l+niE.-E,)} In ^ j 

+ ^ dS, WE,)][l+n(V^)l In ( e -^_ e -^ l+g ,) j } • (G12) 
Changing variables from Ei to u = e~ f5El , the above can be rewritten as 

x 



[Jy \l-u u-yj \u{l-x)J 

+ /•*,( J- + J-kf!!fi^n, (G13) 

Jo \1— u u—y J \ y—ux j J 



110 



where 



y = e- pE « , (G14) 
x ee e"^ . (G15) 



Another change of variable 



z=?P-, (G16) 



1-u 

and the definition of dilogarithmic function, 

Li 2 {z) ee / — — In(l-zu) = £ - » (G17) 
7o u ^ n 2 

yield, 

ST 2 flql 

s ® = ^ 1+n ^l d|k l 



y J VI— x/ VI— x 

In the |q| — > limit, x = y. Hence, immediately, 

ST 2 

S(m th , 0) = — - [l+n(m th )] Li 2 (e"^) , (G19) 

OZ7T J 

and 

Sr°- loop (m th , 0) = ^ Lz 2 (e^ m -) . (G20) 
When the temperature is high, T 3> m th , the identity 

7T 2 

Li 2 (x) = \n(x) ln(l-x) - Li 2 (l-x) , (G21) 

6 



yieldF] 



22 



Similar result was obtained by Parwani pi 
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s twoioo P(mth) Q) = W + 0{mth/THmth/m . (G22 ) 
7b87r 

When the temperature is low, T <C m th , the dilogarithmic function Li 2 (e~^ mth ) = e~ mthf3 , 
and 

Sf °- Ioop (m th , 0) = — e~ m ^ x (1 + 0(e~ m ^)) . (G23) 

1287T 

In the opposite limit where |q| — > oo, the only term that survives in S is the first term 
in the bracket 

l™%) = -3^! <i|k| ln(l-e-«) . (G24) 
In the high temperature limit, this yields 

lta%) = l^(y + 0(>^/r 2 )) , (G25) 



and 



Jim Sr lo ° P (g) = ^(1 + ^K 2 h /T 2 )) . (G26) 



Note that the leading order terms are the same. In the low temperature limit, 



lim Ef " 10015 ^) = -4r- e_mth/3 x (1 + 0(e- m ^)) . (G27) 

<?-°° 1 - 128tt V 2tt 3 v v ;; v ; 
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Numerical integration results shown in Fig. ^ confirms this. A closer look at £/((?) 



with m t h/T = 0.01,0.1 for small values of |k|/T is presented in Fig. p?L and for 



m t h/T = 1.0,2.0,3.0 is given in Fig. |38|. Note that as the mass increases, the self-energy 
becomes flat throughout the momentum range. 
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FIG. 38. Numerical results for £i(E k ) = 7687rE/(£ , fc )/A 2 T 2 for m th /T = 1.0,2.0,3.0. 
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